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Abstract 



In this thesis, conducted in the field of few-body nuclear physics, I studied low energy 
inelastic reactions of electro-weak probes on light nuclei. 

Such interactions hold experimental and theoretical attention in recent years be- 
cause of the role they play in various fields of physics, stretching from checking the 
standard model limits; through nuclear structure and dynamics research; up to mi- 
croscopic interaction in astrophysical phenomena, such as supernova explosion and 
the nucleosynthesis of the elements. 

The reactions considered in this work are: 

(i) Neutrino interaction with A = 3 and ^He nuclei, in energies typical to core- 
collapse supernova. For ^He both the neutral and charged current reactions 
were calculated, whereas only neutral current reactions are calculated for A = 3 
nuclei. 

(ii) Photoabsorption on ^He. This long standing problem has been the attraction 
of theoretical and experimental efforts in the last three decades. 

A complete theoretical description of these reactions holds many obstacles. The most 
significant one arises from the fact that light nuclei are weakly bound and as such 
they have few, if any, excited states. Thus, evaluation of the nuclear dynamics usually 
includes a transition between the nucleus ground state and a resonance in the contin- 
uum, namely fragments of the initial nucleus interacting with each other in what is 
commonly known as final state interaction (FSI). As a result, a theoretical study of 
the reaction should contain both the ground state and the final state interaction. In 
order to avoid a calculation of a continuum final state, which is currently out of reach 
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already for A = A, I used an integral transform with a lorentzian kernel. This novel 
method, the Lorentz integral transform (LIT) method allows a reduction of the full 
scattering problem into a Schrodinger like equation with boundary conditions of a 
bound state. This equation is solved using the recently developed method of effective 
interaction hyperspherical harmonics (EIHH). The combination of these two methods 
has been applied in the literature for electro-magnetic and strong reactions. In my 
thesis, the approach was applied for the first time for weak reactions, and results in 
a percentage level numerical accuracy. 

Photon reactions on nuclear targets are a common experimental method for inves- 
tigating the nuclear structure and dynamics. In electromagnetic processes, the nuclear 
forces manifest themselves also as exchange currents, due to gauge invariance. How- 
ever, one does not have to explicitly calculate their contribution. The conservation of 
the electromagnetic current allows an application of the Siegert theorem, which in low 
energy suggests that these currents are included in the single nucleon charge opera- 
tor. Thus, the scattering operators are model independent. This makes photonuclear 
processes ideal for experimental study of the nuclear force and structure. 

A fundamental example of electromagnetic response which drew continuous in- 
terest in the last three decades, both in theory and in experiment, is the ^He pho- 
todisintegration process. The a-particle is drawing such a great attention because it 
has some typical features of heavier systems (e.g. binding energy per nucleon), which 
make it an important link between the classical few-body systems, viz deuteron, triton 
and ^He, and more complex nuclei. For example in ^He one can study the possible 
emergence of collective phenomena typical of complex nuclei like the giant dipole 
resonance. Furthermore, "^He is the ideal testing ground for microscopic two- and 
three-body forces, which are fitted in the two- and three-body systems. One expects 
that the 3NF is of considerably greater relevance in the four-body system, as the 
number of triplets is bigger, and due to its higher density. 

The microscopic calculation of the ^He photoabsorption process presented here is 
the first evaluation of the cross-section which includes realistic NN potential, Argonne 
1^18, and 3NF, Urbana IX, and full final state interaction, taken into account via the 



LIT method. The resulting cross-section is characterized by a pronounced peak, due 
to a dipole excitation. We have found that the 3NF only slightly lowers the peak (by 
about 6%). This effect is much smaller than the 20% effect on the binding energy, 
and than the influence expected when compared to the ^H/^He(7) processes. The 
effect of the 3NF grows at higher energy transfer to about 35%. 

The comparison to experiments is unclear, mainly due to the poor status of the 
experimental data, in which the differences reach a factor of 2. However, close to 
threshold the theoretical cross section agrees quite well with all experiments. In 
the giant resonance region of the cross-section, where there are several experimental 
evaluations of the cross-sections which in some cases do not agree, the theoretical 
results are in good agreement with the majority of the available data. 

The photoabsorption process can be used to infer properties of the ground state. 
In fact, moments of the total photoabsorption cross-section lead to known sum-rules, 
which are indicative for different characteristics of the ground state, such as the charge 
radius and the nucleon-nucleon distance. In particular, the configuration tetrahedral 
symmetry of ^He is tested, and this symmetry is found to be slightly broken. This 
observation, which can be tested in experiments, suggests that the mean distance 
between identical nucleons is slightly (6%) larger than the mean distance between 
proton and neutron. 

The second branch of the work investigates inelastic neutrino scattering on A = 3 
and ^He nuclei in core-collapse supernova scenario. Inelastic neutrino interaction 
with nuclei can potentially influence several properties of the physics in the super- 
nova. Firstly, they may change the structure of the neutrino signals. Secondly, they 
can deposit energy in the matter behind the stalled shock, thus might change its tem- 
perature and even revive the shock. Finally, they change the chemical composition of 
the star by breaking nuclei to fragments. The high abundance of ^He nuclei both in 
the shocked area and in the outer layers of the star, may result in an effect on these 
processes. 

Lately, we have discovered that a substantial amount of trinuclei can exist in the 
neutrinosphere of the newly born proto-neutron star. This is a result of the previously 
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neglected interaction between nuclei in that area. Thus, neutrino interaction with 
these nuclei may affect the neutrino spectra. 

These facts motivate evaluation of the cross-sections. In this work, the first ab- 
initio calculation of the inelastic reactions is given, with realistic forces: AVI 8 NN 
potential and UIX 3NF. By using the LIT method, all break-up channels are consid- 
ered. The vector meson exchange currents are implicitly calculated using the Siegert 
theorem. The contribution of the (non-conserved) axial meson exchange currents 
to the cross sections has to be calculated explicitly. I cope with this by using ef- 
fective field theory of QCD at low energy, with the nucleons and pions as explicit 
degrees of freedom. The axial currents are calibrated to reproduce the triton half- 
life, within this nuclear Hamiltonian, thus enabling parameter-free prediction of the 
cross-sections. This calculation procedure, known in the literature as EFT*, has been 
tested successfully in the calculation of numerous weak reactions, e.g. Park et. al 
calculated pp-iusion and hep process. The current work is the first application of the 
method for the calculation of complete inelastic processes, not only on-threshold. 

The contribution of the exchange currents to the u — a scattering process is negli- 
gible, due to the "closed shell" character of this nucleus. However, the effect of meson 
exchange currents in the case of neutrino scattering on mass-three nuclei is about 20% 
of the cross-section for low neutrino temperature. As neutrino temperature increases, 
the effect reduces to less than 2% for neutrino temperatures above 5 MeV. 

The error estimation, for both a and the trinuclei, due to the effective theory 
cutoff dependence of the cross-section is less than a percent, for neutrino temperatures 
higher than 2 MeV. This is a strong validation of the calculation. 

The neutrino scattering cross-section is found to be rather sensitive to the prop- 
erties of the force model. Thus, the predicted neutrino scattering cross-sections are 
considered to be of 5% accuracy. 

The reactions calculated in this thesis make an important step towards more 
precise predictions made by nuclear physics. These predictions have the ability to 
check the properties of the nuclear forces on the one hand, and on the other hand 
to be used as accurate and reliable microscopic input for simulations of astrophysical 
phenomena. 
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Introduction 



The physics of the atomic nucleus has drawn a large amount of attention since the 
dawn of quantum mechanics, fascinating physicists in its diversity and complexity. 
Understanding the nuclear world has led to many practical uses, the examples stretch 
from medical application, where radioactive and nuclear magnetic resonance map- 
pings are already standard procedures; through nuclear reactors, in which one makes 
use of the energy released in stimulated nuclear decays as a power source; and finally 
to the destructive power of nuclear weapons. 

The nucleus is a complex entity, built of nucleons held together by the strong 
force. A nucleon is a common name for both neutrons and protons, emphasizing the 
approximated symmetry of their response to the strong force, and their similar mass. 
Determining the properties of the strong force was, and still remains, a scientific 
challenge. It is accepted today that the strong force is the low energy appearance 
of quantum chromodynamics (QCD). However, deriving the nuclear force from the 
fundamental theory is an open question, due to the non-perturbative character of 
QCD in low energy. 

The modern nuclear physics has two main approaches to meet this problem, a 
semi-phenomenological approach and an effective field theory (EFT) approach. In 
the semi-phenomenological approach the potential includes all functional shapes and 
operators allowed by Poincare and isospin symmetry of QCD. The parameters for this 
potential are calibrated by experimental scattering phase shifts. A clear disadvantage 
is that the fit of the parameters is not unique as the theory lacks a clear connection 
to the underlying theory, a fact which leads to a questionable predictive measure. 
However, the approach has many successes. One of its important consequences is 
identifying the need of adding attractive three nucleon forces (3NF) to fit the spectra 
oi A > 3 nuclei. The realistic semi-phenomenological potentials combined with 3NFs 
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successfully reproduce the low lying spectra of A < 12 nuclei. 

The EFT approach uses the chiral symmetry of QCD to describe the force among nu- 
cleons. The approach leads to a consistent perturbative expansion of the interaction, 
with numerical coefficients which contain information about higher energy degrees 
of freedom. These coefficients can be in principal calculated from QCD, but due to 
the aforementioned difficulties are usually calibrated using low energy experiments. 
Among obvious advantages, in this approach the 3NF are predicted and appear at 
higher order of perturbation theory, thus less important - as also found phenomeno- 
logically. However, in order to reach the accuracy achieved in describing the entirety 
of problems as the semi-phenomenological approach one has to expand the poten- 
tials up to next-to-next-to-next-to leading order (N'^LO). Such potentials are still 
not completely available (see however (Epelbaum et al., 2000; Entem and Machleidt, 
2003)). 

It is of great interest to check the success of both approaches in predicting dy- 
namical observables of nuclei. This is not a simple task. The problem arises from the 
need to solve the Schrodinger equation in 3A — 3 dimensions, with substantial inter- 
particle correlations. Without any further approximations, this is an example of the 
few-body quantum problem, which explores the response of A bodies to the forces be- 
tween them. Even today, this problem lacks a general mathematical solution, mainly 
in cases of scattering to continuum states. For bound state problems, the last decade 
has brought immense progress, due to theoretical advances and increasing available 
computer power. A variety of methods solving the nuclear problem microscopically 
(for bound states) have evolved, among which are the No Core Shell Model (Navratil 
and Barrett, 1996), Green Function Monte Carlo (Wiringa et al., 2000), Stochastic 
Variational Method (Suzuki and Varga, 1998) and the Effective Interaction Hyper- 
spherical Harmonics approaches, which have been used for calculations of nuclei in 
the medium range {A > A). 

Contrary to this, calculations of reactions have not met this success. Specifically, 
light nuclei are weakly bound and as such they have few, if any, bound excited states. 
Thus, experiments designed to study nuclear dynamics include a transition between 
the nucleus ground state and the continuum. The latter are essentially fragments 
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interacting with each other in final state interaction (FSI). It is clear that a complete 
theoretical evaluation of the reaction should include an accurate description of both 
the ground state and the final state interaction. An exact calculation of a final state 
in the continuum for all break-up channels is out of reach already for A = 4. In order 
to avoid a calculation of a continuum final state one can use an integral transform 
with a lorentzian kernel. This novel method, the Lorentz integral transform (LIT) 
method (Efros et al., 1994) reduces the full scattering problem to a Schrodinger like 
equation with bound state boundary conditions. This equation can be solved using 
the above mentioned few-body methods, and in this thesis the modern method of 
effective interaction hyperspherical harmonics (EIHH) method (Barnea et al., 2000, 
2001b; Barnea and Novoselsky, 1997) is used for solving the LIT equations. 

The combination of the LIT and EIHH methods allows nuclear theory to describe 
reactions of complex nuclei with A> 3 which were previously unreachable in ab initio 
calculations. The methods have proven successful for the calculation of the total pho- 
toabsorption cross-section (Orlandini, 2004) and inclusive electron-scattering process 
oi A = 3 nuclei (Efros et al., 2004), with realistic nucleon-nucleon potential and 3NF. 
Calculations of photoabsorption and electron scattering processes using semi-realistic 
forces were done for A = 4 — 7 body systems, giving quite realistic results in compar- 
ison with experiments (Bacca et al., 2004, 2002; Efros et al., 1997b, c). 

The focus of this thesis is the accurate modeling and calculation, feasible by using 
the LIT and EIHH methods, of electro-weak interaction with light nuclei, specifi- 
cally neutrino scattering and photoabsorption on nuclei. Electro-weak reactions are 
of central importance not only due to their use in experiments which probe nuclear 
structure and force properties. Their importance stems also from the role they play in 
cosmic and stellar events. This research field, "nuclear astrophysics", evolved in the 
last century from the understanding that in many astrophysical phenomena the typ- 
ical thermodynamic conditions are so extreme that nuclear reactions can be induced 
thermally. By this observation, nuclear physics has turned to be the microscopic 
ingredient used to address astrophysical mysteries. For example, energy released in 
nuclear fusion has been discovered to be the fuel of stars, thus nuclear theory sets the 
main properties of a star, i.e. its composition, evolution and way of death. Another 
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key question addressed in nuclear astrophysics is the origin of the variety of elements 
in nature. Big-bang nucleosynthesis, i.e. the synthesis of light elements in the early 
stages of the universe, provides a hatch to physics in the first three minutes of the 
universe. The rest of the elements are believed to be produced during the evolution of 
a star and mainly in the extreme conditions of its death, in a supernova. In this work, 
neutrino reactions with A = 3 and ^He are calculated. These reactions are of interest 
in the description of the exploding death of a massive star, and in the nucleosynthesis 
within, as will be explained. 

In the next two sections we will elaborate on the specific processes addressed in 
this thesis. 

1.1 Photodisintegration of ^He 

Photonuclear reactions are a common experimental method for investigating the nu- 
clear structure and dynamics. In these processes, a real photon excites a nucleus. 
The perturbative character of the electromagnetic interaction, dictated by the fine 
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structure constant |^ -C 1, makes it possible to seperate the scattering process and 
the nuclear properties affecting the process. The photon, as a result, is ideal for 
experimental study of nuclei. 

A fundamental example of an electromagnetic response which drew a continuous 
interest in the last three decades, both in theory and in experiment (see (Nilsson 
et al., 2005; Shima et al., 2005) and references therein), is the ^He photodisintegration 
process. The a-particle is drawing such a great attention because it has some typical 
features of heavier systems (e.g. binding energy per nucleon), which make it an 
important link between the classical few-body systems, viz deuteron, triton and ^He, 
and more complex heavy nuclei. For example, in ^He one can study the possible 
emergence of collective phenomena typical of complex nuclei. Such a phenomenon 
is the giant dipole resonance, an enhanced cross section in the range 10 — 30MeV 
found in heavy nuclei, due to an E^ dipole excitation. Furthermore, ^He is the ideal 
testing ground for microscopic two- and three-body forces, which are fitted in the 
two- and three-body systems. At present the 3NF is not yet well determined, thus 
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it is essential to search for observables where it plays an important role. Because of 
gauge invariance, in electromagnetic processes nuclear forces manifest themselves also 
as exchange currents, which have turned out to be very important in photonuclear 
reactions and hence 3NF effects might become significant. For the three-nucleon 
systems photonuclear processes have already been studied (Efros et al., 2000; Golak 
et al., 2002, 2005), e.g. in Efros et al. (2000) it was found that the 3NF leads to 
an almost 10% reduction of the electric dipole peak and up to 15% enhancement at 
higher energy. One expects that the 3NF is of considerably greater relevance in the 
four-body system, since it involves six nucleon pairs and four triplets, compared to 
three pairs and just one triplet in the three-nucleon systems. In addition, the higher 
density of ^He might enhance the importance of the 3NF. 

The ^He(7) reaction represents a very challenging theoretical problem due to the 
fact that the nucleus has no bound excited states, so the full four-body continuum 
dynamics and all possible fragmentation have to be considered. It is of no surprise 
that the current theoretical situation of the ^He photodisintegration is not sufficiently 
settled. Calculations with realistic nuclear forces have not yet been carried out, and 
exist only for semi-realistic nucleon-nucleon (NN) potentials. In Refs. (Efros et al., 
1997b; Barnea et al., 2001a; Quaglioni et al., 2004) it has been shown that such models 
lead to pronounced peak cross sections, in rather good agreement with the data of 
Nilsson et al. (2005) and much different from what was calculated earlier (Ellerkmann 
et al, 1996). 

A similar status holds for the experimental research of the process, as most of the 
experimental work has concentrated on the two-body break-up channels ^He(7, n)^He 
and '*He(7,p)^H in the giant resonance region, but a large disagreement still exists in 
the peak. In fact in two very recent (7,n) experiments (Nilsson et al., 2005; Shima 
et al., 2005) one finds differences of a factor of two. 

It is evident that the experimental and theoretical situations are very unsatisfac- 
tory. In this thesis I present an important step forward on the theory side performing 
a calculation of the total photoabsorption cross section a^ of ^He with a realistic 
nuclear force, namely Argonne V18 (AV18) NN potential (Wiringa et al., 1995) and 
the Urbana IX (UIX) 3NF (Pudliner et al., 1997). 
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1.2 Neutrino Scattering on Light Nuclei in Super- 
nova 

Core-collapse supernovae (SN) are giant explosions of massive stars, above 9 solar 
masses, that radiate 99% of their energy in neutrinos. The current theory of core 
collapse supernova holds some open questions regarding the two important phenom- 
ena related to the event, i.e. the explosion mechanism and the synthesis of complex 
nuclei. 

The extreme conditions within the supernova make nuclear reactions important 
microscopic phenomena which govern the equation of state. Therefore, the dynam- 
ics and neutrino signals can be sensitive to the details of neutrino interactions with 
nucleonic matter. In order to analyze the mentioned questions, a better understand- 
ing of the involved microscopic processes is needed. In particular, due to the high 
abundance of a particles in the supernova environment, the inelastic neutrino-^ He 
reaction has drawn attention in recent years. This interest yielded a number of stud- 
ies trying to estimate the cross-section and the role of neutrino-^He reactions in the 
described phenomena (Haxton, 1988; Epstein et al., 1988; Meyer, 1995; Yoshida et al., 
2005, 2006; Ohnishi et al., 2006; Suzuki et al, 2006; Woosley et al., 1990). However 
to date, a full ab-initio calculation that includes a realistic nuclear Hamiltonian is 
still missing. Moreover, the contribution of meson exchange currents (MEC) to this 
particular scattering process was never estimated. 

Core collapse supernovae are believed to be neutrino driven explosions of massive 
stars. As the iron core of the star becomes gravitationally unstable it collapses until 
the nuclear forces halt the collapse and drive an outgoing shock. This shock gradually 
stalls due to energy loss through neutrino radiation and dissociation of the iron nuclei 
into a mixture of a particles and free nucleons. 

At this stage, the proto-neutron star (PNS) cools mainly by emitting neutrinos in 
enormous numbers. These neutrinos are a result of thermal pair production, and thus 
are produced in flavor equilibrium. The characteristic temperatures of the emitted 
neutrinos are about 6 — 10 MeV for i^^^^ (^'a'.t), 5 — 8 MeV for z/e, and 3 — 5 MeV for 
Ue- The difference in temperature originates from the large cross-sections for z/g, z/g 
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electron scattering and charge current reactions, thus u^^r ipiJL,T) decouple deeper 
within the star, where the temperature is higher. Interactions of neutrino with nuclei 
in this location, called neutrinosphere, can change neutrino temperature and spectra. 
Clearly, it is essential to estimate the abundances of nuclei in that region. 

Abundances of nuclei are usually predicted using nuclear statistical equilibrium 
(NSE) models based on binding energies and the quantum numbers of nuclei. How- 
ever, NSE models only treat approximately (or neglect) strong interactions between 
nuclei, and consequently break down near the neutrinosphere. In a recent work 
(O'Connor et al., 2007), we have shown this results in a substantial amount of A = 3 
nuclei near the newly born proto- neutron star (PNS), an effect forehand unnoticed. 

The hot dilute gas above the PNS and below the accretion shock contains up to 
70% ^He nuclei. The high temperature of heavy flavored neutrinos migrating out of 
the PNS leads to a considerable amount of /i and r neutrinos (and anti-neutrinos) 
which carry more than 20 MeV, hence may dissociate the ^He nucleus through inelastic 
neutral current reactions. If these reactions deposit enough energy in the matter 
behind the shock, they can eventually reverse the flow and revive the shock. This 
delayed shock mechanism, originally introduced by Colgate and White (1966), has not 
yet been proven in full hydro-reactive simulations. Haxton (1988) has suggested that 
inelastic neutral reactions of neutrinos with ^He can lead to an enhanced neutrino 
energy deposition. This effect is usually ignored (see however (Ohnishi et al., 2006; 
Woosley et al., 1990)) and was not considered in full hydrodynamic simulations. 

The energy deposition also creates the needed conditions for the r-process, nucle- 
osynthesis by rapid capture of free neutrons, believed to occur in the material ejected 
from the PNS. The break-up of ^He by neutrinos is part of the chain of reactions 
which determines the amount of free neutrons (Fuller and Meyer, 1995) needed for a 
successful r-process. 

A different nucleosynthesis process which is influenced by z/ — a interaction is 
z/-nucleosynthesis in the outer layers of the star. A knock out of a nucleon from a 
^He nucleus in the helium rich layer, creates the seed to light element nucleosynthesis 
in the supernova environment (Woosley et al., 1990). Followed by a fusion of the 
remaining trinucleus with another a particle, this will result in a 7-body nucleus. 
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This process is an important source of ^Li, and of ^^B and ^^F through additional 
a capture reactions. Due to the high dissociation energy of the a, this mechanism 
is sensitive to the high-energy tail of the neutrinos. Thus a correct description of 
the process must contain an exact, energy dependent cross-section for the neutral 
inelastic a — v reaction, which initiates the process. The relatively low temperature 
of the z/g's and i^e's emitted from the core of the star suppress the probability for 
inelastic reactions of these neutrinos with ^He in the supernova scenario. Oscillations 
of the /i and r (anti) neutrinos can yield a secondary source of energetic electron 
neutrinos. The resulting charge current reactions, which have larger cross-sections 
than the neutral cross-sections, would affect the discussed yields (Yoshida et al., 
2006). 

The work in this thesis predicts, to a percentage level accuracy, the cross-sections 
for neutrino interactions with A = 3 and ^He nuclei. This removes the uncertainty, 
which exists today, in the microscopic information needed for a correct evaluation of 
the aforementioned processes. 

1.3 About this thesis 

The work presented in this thesis summarizes the results published in a series of 
papers written during my PhD. studies. 

Two of these papers investigated the photoabsorption process on ^He: 

PI "Photoabsorption on ^He with a Realistic Nuclear Force", D. Gazit, S. Bacca, N. 
Barnea, W. Leidemann, G. Orlandini, Phys. Rev. Lett. 96, 112301 (2006). 

P2 "Photonuclear sum rules and the tetrahedral configuration of ^He", D. Gazit, N. 
Barnea, S. Bacca, W. Leidemann, G. Orlandini, Phys. Rev. C 74, 061001 (R) 
(2006). 

A substantial mass of the work was done on the subject of low-energy neutrino 
scattering on light nuclei, namely triton, ^He and ^He, published in: 

P3 "Neutrino neutral reaction on ^He: Effects of final state interaction and realistic 
NN force", D. Gazit and N. Barnea, Phys. Rev. C 70, 048801 (2004). 
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P4 "Low-Energy Inelastic Neutrino Reactions on ^He", D. Gazit and N. Barnea, 
Phys. Rev. Lett. 98, 192501 (2007). 

P5 "Neutrino breakup of A = 3 nuclei in supernovae", E. O'Connor, D. Gazit, C. J. 
Horowitz, A. Schwenk, N. Barnea, Phys. Rev. C 75, 055803 (2007). 

P6 "Few body Calculation of Neutrino Neutral Inelastic scattering on ^He" , D. Gazit 
and N. Barnea, Nud. Phys. A 790, 356 (2007) . 

P7 "Low energy inelastic neutrino reactions on light nuclei" , D. Gazit and N. Barnea, 

in prep. 

The structure of the thesis is as follows. The first part, chapters 2-5, contains 
a detailed description of the underlying physics and the methods of calculations. 
Chap. 2 explains the quantum formalism of electro-weak probes scattering on nuclei. 
The next chapter uses EFT approach to derive the scattering operator. The LIT 
method is discussed in Chap. 4, followed by a description in Chap. 5 of the nuclear 
wave functions used to solve the LIT equations. 

The second part of the thesis presents the calculated cross-sections and their 
physical application. In Chap. 6 the photoabsorption on ^He process is analyzed, 
based on the work published in PI and P2. Chap. 7 a complete calculation is given 
for neutrino inclusive scattering on A = 3 and ^He nuclei, as was first presented in 
P3-P7. 

The thesis concludes with a summary. 






Scattering of electro-weak 
PROBES ON Nuclei 



The unification of the electro-magnetic and weak forces is considered one of the biggest 
accomplishments of physics in the twentieth century. In the forty years since the 
pioneering work by Weinberg, Salam and Glashow, the standard model has been a 
matter of constant interest and ongoing research. 

The application of this theory to the description of photon and neutrino scattering 
off matter is a key ingredient in experiments checking the limits of the standard model 
or exploring the structure of baryonic matter. These reactions are also the microscopic 
engines of many stellar phenomena, from star evolution to supernova explosion. 

The weakness of the electro-weak interaction, in comparison to the strong force, 
makes a perturbative treatment of scattering processes on nuclei a good approxima- 
tion. Consequently the interaction between electro-weak probes and nuclei targets 
reduces to a current-current type, as shown in the first section of this chapter. 

An important result, presented in this chapter, is that the observables directly re- 
late to the nuclear matrix element. Not only does this simplify calculations, but also 
provides an experimental separation of the electro-weak effects and the strong force 
effects. One cannot overemphasize the importance of the latter fact in the experimen- 
tal research of the strong force, as reflected in the nuclear structure and dynamics. 
From the point of view of a theoretician, a construction of the electromagnetic and 
weak currents within the nucleus is needed. Chap. 3 is dedicated to meeting this 
challenge. 

An additional simplification is achieved by using a multipole decomposition of 
the nuclear currents. This breaks the current into a series of operators, each charac- 
terized by known parity and angular momentum J'^ . For low energy probes, whose 
wavelength q~^ is much larger than the typical size of a nucleus i?, the multipole 
decomposition provides a polynomial dependence in the small parameter qR. In the 
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following discussion, the multipole decomposition technique is rigourously reviewed 
for the electro-weak interaction. 

2.1 Neutrino Scattering 

Neutrino scattering processes on elementary particles can be easily calculated using 
the Feynman diagrams of the standard model. Nuclei, however, are not elementary 
particles. On the contrary, they are built of protons and neutrons, that are built of 
quarks. 

A possible way to cope with this problem, which is common in the literature, is 
to treat the nucleons as point fermions, with form factors indicating their complex 
quark structure. 

The weak coupling constant G = 1.166 x 10^^^ MeV~^ provides a small parameter, 
so the relevant neutrino scattering processes include one boson exchange. The first 
order Feynman diagrams for neutrino scattering off nucleon appear in Walecka (1995). 
The neutrino scattering processes of interest appear in Fig. 2.1. A neutrino with 
momentum fcf interacts with a nucleus of initial momentum Pj^, through the exchange 
of heavy boson with momentum q'^ = {uj,q). The result of the interaction is the 
emergence of a lepton with momentum k't, and nuclear fragments with momentum 
Pr*. Energy momentum conservation implies: g^ = A;f — k't = P't — Pl^. 

The lepton could be either neutrino of the same flavor in the case of neutral 
current reaction, or an electron, muon or tau, for charged current reactions. In the 
low energy regime typical to supernova environment, the possible charged reaction is 
the creation of an electron or a positron. 

Thus, for all neutrino scattering processes, one can write the general form of the 
transition matrix for both charged and neutral currents (with one boson exchange): 

Tfi = (neutrino current)'^ ■ — ^ — |- ■ (nuclear current)'^. (2.1) 

For low energy reactions with respect to the boson masses - M^o ~ 92 GeV and 
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Figure 2.1: Neutrino scattering processes and their kinematics: (a) charged current 
e"*" production; (b) charged current e~ production; (c) neutral current process 

M\Y± ^ 80 GeV - the momentum dependent terms are neghgible. The boson prop- 
agator is constant in this hmit ~ g^^/M"^, and acts as an effective couphng con- 
stant (which explains the wealc scale of weak interaction). This effective vertex, 
demonstrated in Fig. 2.2, leads to the famous current-current Hamiltonian density: 
'Hw ~ J^iJfi U^ is the lepton current, and J^ is the nuclear current). The standard 




c> 




Figure 2.2: Effective vertex for low energy weak interactions. 

model gives additional information regarding the formal structure of the currents. It 
is shown in Appendix A, that the weak neutral current takes the form. 



J. 



(0) 



l-2-sin2M?Jr + ?J,f 



2-sin^ew-J^ 



1 



(2.2) 



To 



|(j; + j;^)-2.sin2^^j; 



em 
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Where J^ and J^ have vector and axial symmetry, respectively. The second equality 
takes advantage of the conserved vector current hypothesis, by which the weak-vector 
current is an isospin rotation of the electro-magnetic current (J!;™ = ^-^JY)- 
A similar calculation for the charged current yields: 

J^^' = YiJ: + J^) (2.3) 

The vector and axial currents include not only a sum over all nucleons. It is well 
known, and will be shown in the next chapter, that in order to keep consistent with 
the chiral symmetry one has to include not only one nucleon currents but also meson 
exchange currents, which hold the nucleons bound in the nucleus. 



2.1.1 Cross section calculation for neutrino scattering pro- 
cesses 

In order to evaluate the probability for a scattering, the starting point is the Golden 
Rule (by Fermi). Taking into account the recoil of the nucleus, the differential cross- 
section takes the form, 



da = 2njde6{e - uj + ^)^s^\imw\^)\'S{E, - E, - e) (2.4) 



Here Ma is the initial nucleus mass, |z) (|/)) is the initial state of the system, and Ei 
(Ef) is the initial (final) energy of the nuclear system. In the previous chapter it was 
shown that, in the limit of low energy-momentum transfer with respect to the heavy 
gauge bosons W^ and Z^ mass, a proper approximation for the weak interaction 
between leptons and baryons is, 

n^ = -^jd^xj^{x)J^{x). (2.5) 

The calculation of the matrix element of the Hamiltonian needed in Eq. (2.4) is divided 
accordingly. Keeping the theory to first order in the weak interaction constant, the 
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lepton current jfj,{x) is composed only of one lepton, which is a point Dirac particle. 
Thus, its matrix element is {f\j^iix)\i) = /^e"*'''^, where for neutrino reaction /^ = 
'u(/ci)7^(l — 75)m(A;2), and for anti-neutrino reaction /^ = 'u(— /c2)7/i(l — l5)v{—ki). 
Consequently, 

{f\nw\^) = ^ j d'xe-''^-^[loJ%{x) -I- Jf^ix)]. (2.6) 

To continue, it is useful to choose the z direction parallel to the momentum transfer 
g, so the plane wave can be expanded as 



^-iq-x 



^zV4vr(2J+ l)jj(gx)r7o(-x), (2.7) 



J 



where jj{x) is the spherical bessel function of order J, and Yjm{x) are the spherical 

harmonics. 

Since the vector / is written in the same spherical orthonormal coordinate system 

e±i,eo = - it is useful to invert the definition of the vector spherical harmonics: 

YjLM = [Yl ® e]i^-J\ to get Yl^cx = Ejm {LmlX\JM) Yjlm- 

For the nuclear current J^^{x), a multipole analysis of the current provides natural 

expansion and some insight in low energy regime. We define the multipole operators, 

viz the Coulomb, electric, magnetic and longitudinal operators: 

CjMiq) = f dxjj{qx)YjMix)Jo{S) (2.8) 

EjMiq) = - [dxVx [jj{qx)YjjM{x)] ■ J{x) (2.9) 

q J 



MjMiq) = J dxjjiqx)YjjM{x) ■ J{x) (2.10) 

LjMiq) = - I dxV[3.j{qx)YjM{x)] ■ J{x) (2.11) 

qJ 

respectively. 

For unoriented and unobserved targets, one has to sum over all possible final 
states, and average over all initial states. For the nuclear sector, this procedure gives 
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the following expression: 



i:^,Y.Y.\m.m(E,^E.-e) = -^^. (2.12) 
E 



X 



1 



-(f- r - /s/s) (7^^, + 7^MJ - 77 (^"^x i').^R^'^E,M, 



+ 



J>0 



Ji is the total angular momentum of the initial nuclear state (which is usually the 
ground state). We use here the nuclear response functions 



'^Ol02 



^6{Ef-E,-e){f\d\\z){f\\d2\\i), (2.13) 



which should be calculated using a model for the nucleus. The^^ sign in the response 
function is a sum over discrete bound states, and an integral over final states which 
are in the continuum. For brevity we use the notation TZq = TZqq- 
For the lepton sector, one sums over all possible helicities of the ejected lepton. In 
the case of neutrino scattering, there is no need in averaging on initial helicities as 
there is only one (anti) neutrino helicity in nature that interacts via the electro-weak 
force. The lepton summation is simply converted to the calculation of traces of Dirac 
spinors, as explained in Appendix B. 

An additional effect, which needs to be taken into account for charged lepton pro- 
duction, is the deformation of the lepton wave function due to the long-range elec- 
tromagnetic interaction with the nucleus. This is done in an approximated way using 
the Fermi function (Fermi, 1934). 

The final result for the differential cross section of a beam of neutrinos scattering 
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on unpolarized and unobserved targets, into solid angle dQ is thus: 

-^^— = — F^''\Zf, kAklx (2.14) 

X I [l - (^, ■ q) {0f ■ q)] J2{n^^ + n,^J Tq-{k- 0f) Yl 2^e7^^,M; + 
I j>i j>i 

+ Y^^]^l-k-Pf + 2{k- q) {Pf ■ q)] 7^|^, + (l + ^. ■ Pf) 7^^, - 2g ■ (k, + Pf) i?e7^e^^._ 



J>0 



''f 



where Zj is the charge of the final nucleus, and Pf = , / ^ is the outgoing lepton 
velocity [rtij is the mass of the lepton). The — sign is for neutrino scattering, and 
the + sign is for anti-neutrino. The superscript (a) = 0, ±1 indicates the process: 

a = Neutral current, for which the G*^"-' = G and the screening function F^^\Zf, kj) = 
1. 

a = ±1 Charge changing current. G^"^^ = G cos 9c ~ 0.974G {9c is the Cabibbo 
angle, mixing the u and d quarks) and the screening function is the Fermi 
function (Fermi, 1934): 

F'^^HZf^kf) = 2(1 + 7o)(2fc;/g^)^(^»-^) I ^[270^+ 1) '|2 ^^' (2-15) 

with a ~ 1/137 is the fine structure constant, Rf is the radius of the final 



nucleus, p = ^aZf/Pj, and 70 = a/1 — («^/)^ (^(a;) is the Gamma function). 

2.2 Relations Between Multipoles and the Siegert 
theorem 

The multipole decomposition is a systematic expansion in the possible angular mo- 
menta exchange in a reaction. This makes it a powerful and intuitive tool to analyze 
reactions, in which the quantum numbers of the final state are given. On a different 
perspective, one can view the multipole decomposition essentially as expansion in 
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powers of the momentum transfer, since jj{p) — (2f+i)\\ ^ot p <^ 1. Thus, it pro- 
vides a tool when investigating low energy reactions with qR < 1 (i? is the radius 
of the nucleus) , for which a finite number of multipoles are sufficient to estimate the 
cross-section to given accuracy. 

The electro-weak reactions discussed in the scope of this work are investigated 
up to pion threshold ~ 140 MeV, thus the multipole decomposition is applicable and 
efficient. 

In this section we will discuss some important properties of the multipole operators 
and relations between multipoles, helpful in calculations. It is useful to understand 
that the longitudinal. Coulomb and electric multipoles are of normal parity, i.e. if 
one denotes the parity of the current by vr then their parity is tt ■ (— )'^. The magnetic 
operator parity is n ■ {—)'^'^^ 

An important relation exists between the electric and the longitudinal operators. 
In order to find it, one uses the following formulas. 



^p[Jj{p)Yjm{p)] = J ^j—j3.J+i{p)Yj J+im{p) + J jj—jjj-i{p)Yj j^im{p) 

Vp X [3j{p)Yjjm{p)] = -iJ^j—^jj+i{p)Yjj+iM{p)+iJ^j—jj.j-i{p)Yjj_iM{p). 



This leads to. 



J + 1 ^ 2J + 1 

Vp X [jjip)YjjMip)] = i^—j—Vp[3j{p)YjM{p)] - -^y J 3j+i{p)Yjj+im{p). 

(2.16) 

Thus, 



EjM = Y — j—LjM - i\l — T — / c?xjj+i(gx)Fjj+iAf (x) ■ J{x). (2.17) 

One notes that since the leading order in Ejm (and Ljm) is 0{{qRY~^), then the 
second term is smaller by a factor of {qRY. Finally, this leads to the following relation 
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for low momentum transfer: 



•^jM ~ y — —LjM- (2.18) 

2.2.1 Siegert Theorem 

An additional relation between multipole operators exists for conserved currents, 
inter alia the electro-magnetic current and the vector part of the weak current. For a 
conserved current V- J = —i[H, JTq], and thus its matrix element (between eigen-states 
of the Hamiltonian) satisfies: V ■ JT" = —iujJq. Let us first look at the longitudinal 
operator: 

LjM{q) = - [ dxV[jj{qx)YjM{x)] ■ J{x) = -- [ dx[jj{qx)YjM{x)] ■ V ■ ^(f) = --CjM{q) 

q J q J q 

(2.19) 
We use Eq. (2.17) to connect the electric multipole to the Coulomb multipole, 



EjM = --\j'^CjM - i\j'^^^ j dx3j+^{qx)Yjj+^M{x) ■ J{x). (2.20) 

These relations are called the Siegert theorem (Siegert, 1937). In low energy, the 
second term in Eq. (2.20) is negligible with respect to the first one. 

The strength of the Siegert theorem is in the connection it makes between currents 
and charge. It implies that for conserved currents in low energy, one can calculate 
all the corrections due to meson exchange currents by calculating the single-nucleon 
Coulomb operator. This not only simplifies calculation, but also keeps the operators 
mo del-independent . 

2.3 Photoabsorption on Nuclei Up To Pion Thresh- 
old 

Photoabsorption is the process in which a real photon is absorbed in a nucleus. The 
process can result in excitation of the nucleus to a bound state, but can also result 
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in a break up of the nucleus, a process called photodisintegration. The latter is the 
case for all photoabsorption processes of "^He, since it has no bound excited states. 

The interaction Hamiltonian of an electro-magnetic process can be found by 
looking at minimal substitution of the momentum in the kinetic energy operator: 
[p — ey4(r))^ = p"^ + e^A^ — e{p, A}. The first two terms are the kinetic energies of the 
nucleus and the photon field, and the last is the interaction term. Hence, it is clear 
that the interaction Hamiltonian is, 

/v 4-7rry f ^ ^ 

Hem = -7^ J d'xA,ix)J''ix). (2.21) 

where a is the fine structure constant, A^ is the photon field, and JT^ is the nuclear 
electro-magnetic current. The photon field in the Coulomb gauge (i.e. V ■ A^(a;)) is 
given by. 



1=^:5: (2..)-/ 



kx^^kx^'"'^ + h.c. 



:2.22l 



1/2 

k A=±l 

Due to the gauge only the transverse polarizations, A = ±1, participate. 

The analogy to the weak Hamiltonian is clear (cf. Eq. 2.5), so the calculations 
are in analogy, but the kinematics is different. 

Let us consider a photon of momentum g and energy uj = \q\, as the photon is 
absorbed. The initial (final) nucleus state is denoted by | i/f), with energies Ei/f 
respectively (by definition uo = Ei — Ef). The photon is absorbed, thus the surviving 
matrix element is (^kX \ at \ O) = ^^/^aa'- 

Due to the possible polarizabilities, only the magnetic and electric enter in the 
calculation. By using the golden rule, we get the final result for the total photoab- 
sorption cross-section: 



a{u) = 47T'auJ2['^EA^)+T^MA^)]- (2-23) 



.Chiral Effective Field 
^Theory 

Nuclear theory is believed to be the low-energy reflection of quantum chromodynani- 
ics (QCD). Though QCD properties are well known, deriving the nuclear forces is not 
possible directly, as QCD is non-perturbative in low energy. 

Two seminal ideas proposed by Weinberg (1979) have pointed the way to a solution 
to the problem. The flrst idea suggests that a Lagrangian which includes all possible 
terms consistent with the symmetries of the underlying theory, would result in the 
most general S-matrix consistent with the theory. The second idea introduces a 
scheme which organizes the terms to dominant and less dominant contributions to 
the S-matrix. This way, the anticipated low energy expansion for QCD becomes 
possible. These ideas were used in the last 15 years to construct a nuclear Lagrangian 
based on an effective fleld theory for low-energy QCD, named chiral perturbation 
theory (xPT) (see for example (Ordofiez et al., 1994; van Kolck, 1994; Epelbaum 
et al., 2000; Entem and Machleidt, 2003)). 

In this chapter, the application of effective fleld theory to QCD is presented, and 
used to get a description of the currents in a nucleus. In turn, these currents interact 
with electro-weak probes, determining their coupling with matter. The effective fleld 
theory is constructed to describe nuclear electro-weak processes of typical energy of 
up to 100 MeV. Thus the relevant degrees of freedom are pions and nucleons. 

3.1 QCD symmetries at low energy 

QCD, as a fleld theory, describes the interaction among quarks and gluons. Its un- 
derlying gauge symmetry is the non-abelian color SU{3) Lie algebra. The masses of 
the up and down quarks, which are the building blocks of nucleons, are small (Yao 
et al., 2006): rriu = 2±1 MeV and ma = 5 ±2 MeV. Thus, these masses are negligible 
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with respect to the QCD mass scale (~ 1 GeV). In the context of nuclear physics, it 
is natural to treat QCD in the limit of zero quark masses, with restriction to the up 
and down quarks. 

Following these restrictions, one writes the quark field composed of two compo- 

nents: I . The massless quark field enters in the QCD Lagrangian as, 

\d J 

Cqcd = ^Ql.^'Q- (3.1) 

The quark field has a left and right symmetry. We define right and left handed 
components of the quark field: q^^L = |(l i Ts)?; respectively. Using the properties 
of the 7 matrices (see Apendix B) , the Lagrangian can be separated, 

'C-lcD = nRlfJ^^R + iqil^D'^QL. (3.2) 

One may check that this Lagrangian is invariant under the global transformations: 



j..exp(-z0«.^)(^ 



qn=\ Ij^^^P(-^®''-i)( 11^ 1 ^3.3) 



and, 



where r are pauli matrices generating the flavor symmetry between up and down. As 
right and left quarks do not mix, this implies the chiral symmetry SU{2)r x SU{2)l 
of QCD. One can rewrite this to a symmetry on the original quark field: 

^ ) ^exp[-zi(l±75)e-^] [ ^ |. (3.5) 
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The chiral symmetry can be transformed hnearly to a vector symmetry, 

«)„..p[_,e^. !-](«), ,3.) 

and axial symmetry, 

Nother's theorem suggests the existence of conserved vector and axial currents: V^ = 
qi^\q and A^ = g7^75|q'. 

The problem is that this symmetry entails a degeneracy of ground states. A vector 
degeneracy is found, the isospin symmetry, which couples for example the proton and 
neutron to a single entity, the nucleon. This is not the case for the axial symmetry, 
as multiplets of opposite parities are not observed in nature. One concludes that in 
nature SU{2)r x SU{2)l is broken to SU{2)v 

This is an example of spontaneous breaking of symmetry, i.e. a symmetry of the 
Lagrangian which is not realized in the ground state. A consequent of thus is the 
existence of Goldstone-Nambu bosons, massless particles that posses the quantum 
numbers of the broken generators. In the case of the chiral symmetry, the Goldstone 
bosons should be pseudoscalar. The pions, pseudoscalar isospin triplet, are considered 
to be the Goldstone bosons of the chiral symmetry. As a result of the non-vanishing 
up and down small masses, the pions obtain mass, but they are light ( ^ ) ~ 0.02. 

The appearance of Goldstone-Nambu bosons is a result of the Goldstone theorem 
(Goldstone et al., 1962). This theorem gives an additional information regarding the 
structure of the transformations. The pion field is defined by, 

^^exp(.|^T ,3.8) 
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Goldstone theorem leads to the transformation rule for U = C,"^ 

U = exp ( ^^ ) ^ LUR\ with L e SU{2)l and R E SU{2)r (3.9) 

Thus, U transforms as the (2, 2) representation of SU{2)r x SU{2)l- This is a non- 
linear realization of the symmetry, as WU = 1 and det U = 1. 

This implies the following transformation rules on the nucleon and pion fields: 

N^hN (3.10) 

e ^ Lih) = hiR^ (3.11) 

with h G S'f/(2)v, local matrix. 

With these transformation laws, it is possible to build Lorentz invariant terms to 
include in the Lagrangian of the effective theory, 

i^eff = Y.CACA (3.12) 

A 

This is an effective theory, thus renormalizability is inserted via the low energy coef- 
ficients. This sum however is infinite. Hence, a logic in deciding the order of a term 
in the expansion is needed. 

3.2 Power Counting of Chiral Perturbation The- 
ory 

EFT of low energy QCD is called Chiral Perturbation Theory (xPT). The theory 
has a natural cutoff scale, the chiral symmetry breaking scale, which is of the order of 
the nucleon mass A ~ 1 GeV. The expansion relies on the small momentum typical 
to the needed process Q, thus the expansion parameter is (Q/A). 

The contribution of an interaction is determined by its power: (Q/A) . Let a 
scattering process be described by a diagram. The diagram is built of C separately 
connected pieces, A nucleons, L loops. The diagram includes several vertices, each 
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involves Ui nucleon fields, di pion masses or derivatives. The order v is given by 
(Weinberg, 1996): 

z/= -2 + 2A-2C + 2L + ^Ai (3.13) 

i 

with Ai = di + ^ - 2. 

One has to note that v is bounded from below, thus enabling systematic expansion. 
This counting holds not only for diagrams, but also for Lagrangian terms. The 
counting for a Lagrangian term is similar to a vertex: i/£ = (i + |. 

The power counting formalism introduced here is applicable when the choice of 
cutoff A is relevant for the energy-regime of the process (see discussions by Nogga 
et al. (2005); Epelbaum and Meissner (2006)). 

3.3 Effective Lagrangian 

As aforementioned, the chiral effective field theory takes into account pions and nu- 
cleons. Not surprisingly, the Lagrangian includes three parts, a pure pion, a pure 
nucleon and pion-nucleon interaction Lagrangian. In this section I will review the 
different parts of the Lagrangian, and emphasize those parts needed for deriving the 
currents within the nucleus. In general, Lagrangians to next-to-leading order (NLO) 
are considered, built using the formalism introduced in Coleman et al. (1969); Callan 
et al. (1969) and following the discussion in Ananyan et al. (2002). 

3.3.1 Pion Lagrangian 

The Lagrangian has to be a Lorentz invariant scalar term. Thus, it should include 
even number of derivatives/pion masses. Since the expansion used here is up to NLO, 
the pion Lagrangian term is expanded to i^ = 2, in the power counting introduced 
previously. Thus, 

42J = ^Tr [d^\]d^\]^ + m^ (U + Ut)] (3.14) 



f2 

I 

/^ ~ 92.4 MeV is the pion decay constant, calibrated from this process at low pion 
four momentum. The pion mass is taken as ??2^ = 139 MeV. 
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3.3.2 Pion— Nucleon Interaction Lagrangian 

The pion-nucleon interaction, to lowest order, arises from the nucleon free Dirac 
Lagrangian. The reason is that the Dirac Lagrangian iN {^^y ^d^ ~ Mq) N is not invariant 
under the transformation rules of Eq. (3.10). This is not surprising, as the Yang-Mills 
theory demands adding gauge term to the derivative. The additional term has to be 
of chiral order 1, equivalent to the derivative. It can be easily checked that this term 
is: 

v^^-'-{^d,i'^ + i'^d,e>- (3.15) 

The chiral covariant derivative is: V^ = d^ + iVf^. The notation f^ indicates the 
vector-like parity of this term. Trivially, a term with axial-like parity of the same 
chiral order is, 

a^^-'-{^d,e-^^d,^). (3.16) 

This leads to the leading order pion-nucleon Lagrangian, 

-^S = N {i-f,V^ + gAl^l^a, - Mo} N (3.17) 

The low energy constant (LEG) introduced here is the axial coupling constant g^ = 
1.2670 ± 0.0035 (Yao et al., 2006). Mq ^ 890 MeV is the nucleon bare mass. 

Expanding Eq. (3.17) to leading order in the pion field, one gets the leading 
interaction Lagrangian: 






^St = -—e'^'^n'^Ni^n'XN + ^N{^n^)j, (3.18) 



Using this, one constructs the Feynman diagrams related to this interaction La- 
grangian. 

Constructing the next-to-leading order contribution to the pion-nucleon interac- 
tion is easier, as the building blocks are a^ and v^, and U. It is easily shown that 
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there are two terms at this order: 

hCf^ = ^NNTT{^,^^) (3.19) 

52C% = i^N[a,, a,]a^''N, (3.20) 

Mm 

Where a^'^ = |[7^,7'^]. The LECs here are cahbrated using ir — N scattering data 
cg = -3.66 ± 0.08, and C4 = 2.11 ± 0.08 (Bernard et al., 1995). 

In view of the need of Feynman diagrams, one gets to this order: 

AS = ^^T . {d,n X d^n)a^^N + -^NNd.i.'^dV (3.21) 

J TV J TV 

3.3.3 Nucleon Contact terms in the Lagrangian 

Nucleon contact terms represent the short-range nuclear force, usually connected to 
the exchange of heavy mesons, such as the u;(782). This mass is too close to the QCD 
scale, so it cannot be described as a multi-pion interaction using xPT. The contact 
terms are important for the purpose of renormalization and dimensional regularization 
of pion-exchange. In case of infinities, the contact terms act as counter-terms. 

The contact terms include two nucleons, without meson fields. Being the represen- 
tative of higher energy modes, these terms are not constrained by chiral symmetry, 
and required to fulfill parity and time- reversal symmetries (Ordonez et al., 1994). 
Parity suggests even numbers of derivatives. The LO has no derivatives, and the 
NLO consist of two derivatives/pion masses. 

Contact terms appear also when pions, electro-weak leptons, photons or other 
probes couple to pairs of nucleons. This gives rise to a contribution of few nucleon 
terms to the nuclear force (Grdestig and Phillips, 2006), and two nucleon exchange 
currents. The leading effects are represented by a term in the chiral Lagrangian (van 
Kolck, 1994), 

£4 = -2/^1 (iV7^75«M^) [NN) . (3.22) 

In general, the LECs in the contact terms should be calibrated by nuclear matter 
experiments. In the current case the short-range LEC is calibrated using the triton 
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half life. This process is discussed and applied in Sec. 7.2. 

3.4 Currents in the Nucleus 

The axial and vector symmetries of the Lagrangian imply conserved Nother currents. 

These Currents are, 

dC 

^'' ^ ~d{d,e'^{x)) ^^-^^^ 

where the superscript a indicates isospin components. e"(x) are local, infinitesimal 
parameters of the symmetry transformation. 

A vector transformation has to be invariant to parity. The definition of parity, P : 
L ^ R, means that for a vector transformation with infinitesimal group parameters 

L = R = expli^] . (3.24) 

The h{x) transformation field can be found using Eq. (3.10), to get h{x) = L = R. 
For an axial transformation one gets 

L = expU^Y R = L^ = expl-t^Y (3.25) 

In this case h{x) should be solved from its definition in Eq. (3.10). Putting h{x) = 
exp [i-^], one can find. 



7 = — 2(/3 X vr) arctan 



tJ3jii Tz tjciri fj J. 



1 — (/9 ■ 7r)tan | tan i^ 



-(/?X7r)tan-- + 0(/?2) (3.26) 



2 """" 2/^ 

Using this, one can derive the currents corresponding to the Lagrangians of the 
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previous section. The vector currents: 

V£)"^ = -iyTr {r" [Ud^'U^ + U^d^U] } = (tt x 9^7f)"sinc2 ( ^ j (3.27) 



3\ 



= (tt X d^^n)" + 0{tt 
VS"^ = ^Nn^r'^^^ + ^r^m + IgANl'-lrA^r-e " ^^r^^ = (3-28) 

b c 

= NY—N[cos{^){6''^ - vr'vr^) + vr^vr''] + ^^e'^''^7r^in(^)iV7''75— iV = 

■^ In ^ 

VS"^" = -^^^"''[^^^'^^ - ^^"^^ <^ - ^A^A^TrieV'^e - er•^e^ a,] (3.29) 
= ^^^'^^ [^ X (S.vf X f )r TV 

where sinc(x) = ^i^^^. The axial currents: 

^(^sj'^A* = -z^Tr {r*^ \Ud^U^ - U^d^U\ } (3.30) 



-/7r<9;,vr 



6 



sine 



— 1(0 — vrvrj+vrTT 



-A9X + 0( 



TX 



•^S"'^ = ^iV7'^[er'^e^ + ev-^eiiv + ^^AiV7^5[er'^e^ - ev-^eiiv = (3.31) 

c fe 

= -e'^*^7r''sin(^)A^7'^— iV - ^AiV7'^75— iV[cos(^)(5"* - tt^tt*) + tt'^tt*] 

^?^^ = -^iVa'-[eV"e + ^r^e, a,]N - ^iViVTr[etr'^e + ^r^^, s] (3.32) 

Af^^ = -DiNYl5T''NNN (3.33) 

The notation used in the upper index denotes: (z/^) is the chiral order, a = 0, ±1 is the 
isospin projection, and /x = 0, .., 3 is the Lorentz vector projection. The expansion is 
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used to construct the Feynman diagrams for the currents (see Appendix D). 

3.5 Nuclear Currents Operators 

Using the equations of the two previous sections, one constructs the Feynman rules. 
The rules derived from the Lagrangian are used to construct nuclear forces form xPT. 
Let us assume a case, in which the interaction Lagrangian is 

L = -gS^{x)J^{x), (3.34) 

where S^ is an outer source current, and J^ is the nuclear current. This Lagrangian 
is relevant for low-energy electro-weak interaction (cf. Eqs. (2.5,2.21)). The Feynman 
rules derived from the Nother of the Lagrangian symmetry, are used to construct the 
nuclear currents. 

An extensive description of the Feynman rules and the derivation of the nuclear 
currents appear in Appendix D. In view of the typical momentum of the processes 
in discussion, one can use a soft probe approximation and write the currents to first 
order in the momentum transfer. It is also clear that any recoil of nucleons is smaller, 
and thus the use of heavy baryon expansion is justified, expanding the nucleon spinors 
to leading order in -^J— (see Appendix C) 

The leading nuclear current operators are single nucleon operators, i.e. describing 
a coupling to only one nucleon within the nucleus, as if it was a free particle. The next 
contribution involves two nucleons, in which case the probe interacts via a contact 
term of two nucleons or a nucleon-pion vertex. 

The resulting matrix elements are expressed in momentum space, and as they 
originate in EFT approach, they are valid up to a cutoff momentum A. The cutoff 
momentum is chosen such that A ^ Q, and lower than QCD scale. A proper choice 
is thus A = 400 - 800 MeV, as in (Park et al, 2003; Ordoiiez et al, 1994; Grdestig 
and Phillips, 2006). The renormalization group theory states that an observable does 
not depend on the cutoff momentum. Any cutoff dependence implies excitation of 
higher energy degrees of freedom. Thus, the cutoff dependence is a consistency check 
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for the theory. In the Lagrangians of the previous section, the only LEG with a 
connection to higher energy modes is dr, thus one should expect a rather steep cutoff 
dependence for this LEG, and its cutoff dependence is chosen in such manner that 
the triton half-life reproduces its experimental value. 

The nuclear wave functions are written in configuration space, thus one should 
use a proper transform to move the operators to this basis. Fourier transform is of 
no help as it includes high momentum. Following Park et al. (2003), let us define the 
cutoff Fourier transform of an n-body current as 

Jl2...n{^l,---,Xn]q) = (3.35) 



n/^^^'"'^A(^i 



{2nf 6^''\q + h + k2 + --- + K)Ji2...nih, ...,K;q). 



The cutoff function, S\(k) = exp(— ^), does not influence one body currents, due 
to the momentum conserving 6 function. For Q ^ A the cutoff function Sx{k) is 
unity, and it approaches zero for Q ^ A. This choice of cutoff function is convenient 
both for numerical and analytical calculations. One should note that this cutoff 
function inserts an additional j^ dependence to the currents. This is, however, a 
minor, negligible effect. 

For low energy reactions, where meson exchange currents (MEG) are a perturba- 
tive effect, one can use the long wavelength approximation q = 0. Thus, the Fourier 
transform of a MEG depends only on the relative coordinate Xi2 = xi — X2, and 

^^^e-'''-^^^Sl{e)AUki = -k, k2 = k). (3.36) 

As the Fourier transform acts on pion propagator, the needed functions are, 

4'^(-1- [l^y'-'Sl{k% (3-37) 



, (27r)3 

[27T)'' ^2 _^ ^2 
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yUr) = -|:l/Xo(0, (3.39) 

yUr) ^ -\r-^-^yUr). (3.40) 

m^ or r or 

One reproduces the usual Yukawa functions in the hmit A — i> oo. 

This process leads to operators acting on one nucleon in the case of single nucleon 
operator, or on a pair of nucleons in the case of MEC operator. To move A nucleons 
(or ^ ~ ' pairs of nucleons) to their coordinates in configuration space one uses first 
quantization. For the single nucleon operators, 



i(^)^(x) - 




-ri), 


and for the MEC opeartors. 






J^^^'^^{x) = 


i<j 


-r,j 



(3.41) 



(3.42) 

This leads to the final form of the operators. The single nucleon operators, which 
are often called the impulse approximation (lA) are the same as those derived from 
standard nuclear physics approach: 

^T'^-..(^,*'"(^-.1)-|{4,{4.*'='(^-^.)}} (3.46) 

\2Mf, \2Mf,' ^ •'// 4Af„\2MA,' * 
Due to the Siegert theorem (Sec. 2.2.1) for low energy reactions, the conservation of 
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vector currents implies that neither the vector current nor the vector MEC should be 
calculated explicitly, as they are included implicitly in the calculation of the single 
nucleon charge operator. 

The axial charge MEC operator is: 

•A%\n,) = -^(f(^) X f(^))'^(a(^) + a(^)) ■ %y,\(r.,), (3.47) 



and the axial MEC, 



-C3(Or - 0^'^)(^5f (rl,) - ^m^yo\)- 

-(C4 + -)&^''mly^^{rij)- 

- (C4 + ^)0^"(^5f (rl,) - ^m^yo\)} + 

+ :^(0^» + 0^»)5f(r-;,) (3.48) 



where 



d» ^ -!^(r(^) X T^'^riP.a^'^ . h2 + P2a^'^ ■ ris) (3.49) 

^0'" = {r\2ri2 - y) C?r, (3.51) 

and © = x,+,-. 

It is important to notice that the contact parts of the current, i.e. terms proportional 

to S{r), do not vanish only for pairs in a relative L = state. This implies S + T = 1, 

meaning either S = or T = 0. Consequently, when O^ acting on these states, it 

vanishes. 

An additional simplification is achieved by noting that O^ and Oq are connected. 
The fermion statistics implies S^S'^S'" = — 1, where S^ exchanges isospin of the two 
particles, 2°" exchanges spin, and S^ exchanges locations. A trivial identity for Pauli 
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matrices: a^ x a^ = i{(r^ - (r'^)E'^ . Thus, O^ = -OqE'^, and as E''6{r) = 1 the 
connection is clear (Park et al., 2003). 

Thus, all the coefficients of S^^\r) are combined into a single dimensionless LEC^ 
dr = — —Di + ^ + |(c4 + j). The 0(A~^) appearing in Sj^ {■f) has minor effect, 
which cancels when renormalizing dr. The renormalization dr{A) for a range of cutoff 
momentum A = 400 — 800 MeV, is calibrated to reproduce the triton half-life. 

Bearing this in mind, the axial MEC are 

42"(^^-) = 77^^>^"^f a) - ?^^pVA(ri2)- (3.52) 



9Aml 

■2M/2 

9Aml 






^^{o^ + o^^) + \{h+^^oi^ 



1 
csl^e' +Te' ) - (c4 + -)T^' 



y2AVij) 



2M/2 

3.6 Standard Nuclear Physics and EFT 

This chapter has shown the vast advantage in applying EFT approach for the nuclear 
problem. The theory provides a perturbative approach with its symmetries completely 
under control. This fact gives a predictive measure to the approach. 

These advantages are missing in the standard nuclear physics approach (SNPA), 
which is phenomenological. However, SNPA approach has many successes, using 
Lagrangians that include Nucleon-Nucleon potentials and three nucleon forces one 
can describe the low lying spectrum of nuclei with A < 12. SNPA Lagrangians are 
used to calculate electro-weak and strong nuclear reactions, to good accuracy. The 
Lagrangians are built to reproduce hundreds nucleon-nucleon scattering processes 
with energies up to 500 MeV, with x^ per datum close to unity. For example, the 
binding energies of light nuclei are reproduced using these forces to accuracy of tens 
of keV. 

In order to reach this level of accuracy with EFT based Lagrangians, one has to 
expand to at least next-to-next-to-next-to-leading order, which includes more than 



-'^ Different variations of EFT (e.g. (Park et al., 2003)) lead to different descriptions of the short 
range forces. These, however, renormalize to the same LEG. 
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20 parameters. This was done only lately, by Epelbaum et al. (2000); Entem and 
Machleidt (2003), and still has not arrived the level of x^ achieved in SNPA. 

Lately, Rho et al. (Rho, 2006; Park et al., 2003) tried to combine the advantages 
of both approaches. The idea relies on the fact that albeit phenomenological, the 
SNPA approach holds correct long range behavior. Thus, one may use the nuclear 
wave functions predicted by SNPA, and calculate reactions using the nuclear currents 
derived in EFT. The typical accuracy for reactions is in the percentage level, hence it 
is sufficient to expand the forces up to next-to-leading order, with only one unknown 
LEG needed for renormalization. 

This hybrid approach, named EFT*, is of course not consistent, as the currents 
and forces are not derived from the same Lagrangian. Nevertheless, it has been 
used to calculate reactions oi A = 2,3,4 nuclei (Rho, 2006; Park et al., 2003). The 
minor response of the calculated cross-sections to changes in the cutoff momentum 
is an evidence for the power of the method. This method will be used in this work 
(Chap. 7) for the calculation of neutrino scattering on yl = 3 and ^He nuclei. 
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The response of a system to excitations is a key ingredient in understanding its 
properties: stability, energy levels, structure, etc. 

For an inclusive reaction with energy transfer u;, described by an interference 
between two operators Oi{q) and 02(g), the response function is defined as: 

^o.o.(^,?") = £'^(^/-^.-^)(^|Ol(g1l/>(/|02(g1K>, (4.1) 

where \i/ f) and Ei/j denote the initial and final wave functions and energies, respec- 
tively. In order to calculate this observable, one has to consider every possible final 
state, including those which are in the continuum. For the light nuclei, which have 
few bound excited states, if any, the continuum must be taken into account if one 
wants to consider uj ^ 0. This makes an explicit calculation of the cross-section for 
such nuclei a difficult task, not only as one has to keep track of all break-up channels, 
but also since a general continuum state wave function is out of reach already for 
A = 4. 

However, the response function is an integral quantity, as it depends only on pro- 
jections of the initial state. This has been the motivation for using integral transforms 
on the response functions, by which, as shown in the next Section, one can avoid an 
explicit calculation of final state. 

One of these transforms is Lorentz Integral Transform (LIT) method. This method 
has been used in the calculation of various strong and electro-magnetic nuclear re- 
actions (see for example (Efros et al., 1997b, c, a, 2000; Barnea et al., 2001a; Bacca 
et al., 2002, 2004)). In this thesis the LIT method is applied for weak reactions on 
nuclei (Gazit and Barnea, 2004, 2007b, a, c). Though used here for inclusive reactions, 
it has been applied for exclusive reactions as well (Quaglioni et al., 2004). 
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4.1 Integral Transformations 

Using integral transform transforms the response function to a continuous function, 
which depends on a parameter a, using an analytic kernel function K{uj,a). The 
integral transform is a functional of the response function, defined as: 



I[n{uj,q)]{a,q) = dcuTZiu , q) K {u , a) (4.2) 

Jo- 



The integration starts from 0~ to include every possible scattering, including elastic 
channels. Finite transform and its existence depends on the high energy (high uj) 
behavior of both the response function and the kernel, and should be checked for 
each case. 

The advantage in using an integral transform is easily viewed when inserting the 
definition of the response function, Eq. (4.1), in Eq. (4.2). The integration over the 
energy conserving 6 function is trivial. 



A<r,q) = ^{z\dMf){fmmKiEf - E,, 



cr 



As K{Ef — Ei, a) is a number, one can insert it to the matrix element. Using then the 
fact that the states |/) are eigenstates of the Hamiltonian, Ti, namely 7i|/) = Ef\f'^, 
yields the expression 

T(a,g-) = f](2|6I(g-)K(7i;-i?„a)|/>(/|62(g1K>- 

One notes that in this expression the final state appears only through the projec- 
tion operator |/)(/|. The hermiticity of the Hamiltonian implies closure:^]; l/)(/| = 
1. Thus, 

X{a,q) = {z\dl{q)K{n-E,,a)d2m^). (4.3) 

There is no doubt that the use of integral transforms simplifies the calculation, 
as it depends only on the ground state properties. However, this is merely a disguise 
as the kernel operator depends explicitly on the Hamiltonian Ti.. The latter acts on 
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02(g) |i), which is not an eigenstate, hence depends on all basis states. An additional 
unresolved question is the inversion of the transform. 

Obviously, an intelligent choice of the kernel can solve both these problems. Nev- 
ertheless, the solution can never be perfect. A choice of a broad kernel smears out 
narrower structures, which when inverting the transform results in numerical insta- 
bility. Choosing narrow kernel leads back to the solution of the continuum state 
problem. 

In the literature, different kernels have been proposed. Among which is the Laplace 
kernel K{u}^ a) = e""'^ . The calculation of the transform with this kernel is relatively 
easy, as it can be presented as propagation in imaginary time. This fact has made 
the Laplace kernel appropriate for Quantum Monte Carlo calculation (e.g. (Carlson 
and Schiavilla, 1994)). The disadvantage is the fact that an inversion of Laplace 
transformation demands an analytic continuation of the integral transform, which is 
usually unavailable. 

In view of the needed narrowness of an optimal kernel, it is valuable to use a 
kernel which depends on an additional parameter aw, with the following properties: 

• K{uj, (7\ (j^j^) = K{uj — a; aw)- 

• lim^j^^o K{uj - a; aw) = S{uj - a). 

The advantages are clear, aw is chosen in view of the anticipated width of the reaction 
cross-section. Photodissociation and neutrino scattering reactions are dominated by 
giant resonance excitations, with a typical width of 10 — 30 MeV. Thus a reasonable 
choice is aw ~ 10 MeV. 

Possible examples for kernels which satisfy the aforementioned properties are the 
lorentzian and the gaussian. The Lorentzian turns to be simple to calculate, hence 
attractive. 

4.2 Integral Transform with a Lorentzian Kernel 

A Lorentzian kernel is defined as K(uj,a) = -, ^, , ;j , where a = —an + iaj is a 

complex continuous parameter, with a real part indicating the transform variable and 
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its imaginary part indicates the Lorentzian width. The Lorentz Integral Transform 
(LIT) is defined as, 

C{a^q)= Hdu ^ '^^''f , . (4.4) 

Jo- iuj-aR)+aj 

Putting the LIT kernel definition in Eq. (4.3) yields: 

£{a,q) = miq) ^ I OMl^), (4.5) 

n — Lo + a* n — Lq + a 

which has the form of an internal product: 

C{a,q) = {^^\^2) (4.6) 

with, 

|v&,)^-^^ 0.(g)H>, for z = 1,2. (4.7) 

This important result implies that the evaluation of the LIT requires the solution 
of the following equation: 

{n-Eo + a)\^iia)) = di\z). (4.8) 

These equations have several appealing properties: 

Uniqueness and existence of the solution The hermiticity of the Hamiltonian 
implies real energies, hence Ef — Eq — a ^ ioi aj ^ 0. Thus the homogenous 
equation has only the trivial solution. As a result, there could only be a unique 
solution dictated by the source term in the right-hand-side of the equation, and 
the boundary conditions are dictated by it. The source term has a finite norm, 
i.e. the boundary conditions of a bound state problem. This means that the 
inner product of Eq. (4.6) is finite and thus the transform exists, and unique. 

LIT is Schrodinger— like problem In view of the bound state boundary condi- 
tions, Eq. (4.8) is a Schrodinger equation with a source term, thus can be 
solved using the variety of methods available for such problems. One of these 
methods is the Effective Interaction in the Hyperspherical Harmonics (EIHH) 
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discussed in the following chapter. 
LIT is a Green— function— like matrix element It is easy to rewrite Eq. (4.5) as, 

C{a,q) = -Im ! {^\^l{q) ^^-^ OMli)] . (4.9) 

Thus, a matrix element of a green-function. Several methods exist to calcu- 
late such a matrix-element, one numerically effective method is the Lanczos 
algorithm, described in the following subsection. 

The Lorentz integral transform clearly simplifies calculation of response functions. 
The method avoids the difficulty of calculating continuum wave functions. Instead 
a bound-state like equation has to be solved, enabling the theoretical evaluation of 
previously unreachable reactions. 

4.2.1 Calculation of LIT using the Lanczos Algorithm 

In this subsection an efficient algorithm for the calculation of the LIT presentation 
of Eq. (4.9) is reviewed, based on the Lanczos Algorithem (Lanczos, 1950). First, let 
us focus on the Oi{q) = 02(g) = O case. 

The Lanczos basis {\ipn);n > 0} is an orthonormal set of states starting with a 
known state \^o}- The rest of the basis can be defined recursively. 



'M+i 



\(Pn+l) = T-Cl^Pn) - an\(Pn) " bn\(Pn^l), for U > 1, (4.10) 



With the condition bo = 0. From Eq. (4.10) it is trivial that the Lanczos coefficients 
a„, and bn are given by: 

an = {'fn\'H\^n) and 6„ = {(Pn\'H\'^n~l) ■ (4.11) 

Choosing the Lanczos pivot as 

Iv^o) = ' ^ , (4.12) 



4^^^^ 
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allows to write the LIT as 

C{a) = -(z|6t6|z)lm{xoo(^o - ^)}, (4.13) 

Where 

Xno{z) = {ipn\:^, IV^O)- (4.14) 

Ti — Z 

In order to calculate this , one uses the identity {7i — z){7i — z^ = I, which leads 
to linear set of equations, 

E(^-^)»^n^"0(^)='^-0• (4.15) 



Applying Kramers rule, one can express Xno{z) as continuous fractions, in particular 

xooiz) = ^ i p . (4.16) 

{z - ao) - ^ '-^ 



(z-a2)- 



(2-«3)- 



In typical examples a few hundreds of Lanczos coefficients are sufficient to reach 
convergence in the LIT. 

After defining Xno{z), a formulation of a LIT calculation with Oi{q) ^ 02{q) is 

available. In this formulation, the Lanczos pivot defined as (/9o) = , ^. One 

y'(i|6|62|i) 

inserts an identity operator J2i \^i){^i\ to get: 



C{a,q) = ^ -lm\2{i\0liq)\ipn)xno{z). (4.17) 

4.3 Inversion of the LIT 

The success of an integral transform highly relies on the existence of a stable inversion. 
An analytic inversion for the Lorentz transform does not exist. 

One can observe the problem in inverting the LIT in the following reason. Let us 
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perform a Fourier transform on Eq. (4.4). The r.h.s of that equation is a convolution 
of forentzian and the response function. Bearing in mind that the fourier transform 
of a forentzian is exponential decay, one gets 

:F[C{a,q)m = -e-'^^l'=l^[7^(c.,g-)](fc). (4.18) 

o-i 

As a result, high energy (short wavelengths) are suppressed when performing the LIT. 
Clearly, this fact creates instability in the inversion, which enhances short wavelength 
oscillations. 

However, a hope lies in the fact that mathematically the inversion is unique. 
Moreover, there are several numerical approaches for the inversion in the literature, 
which have provided stable and successful results. The method used in this scope 
is expanding the wanted inversion (i.e. the nuclear response function) in a set of 
functions /„(co'), 

N 

n{u,q) = Y,^n{q)fniuj). (4.19) 

71=1 

Since the transform is linear, it is obvious that the LIT of Eq. (4.19) C{a,q) = 
J2n=i Cn{(l)'C.n{cr), whcrc £„,(a) is the known LIT of {fn{uj)}. One obtains the coeffi- 
cients Cn by solving the minimal squares problem, on a grid {crjj, i = 1, . . . , K}: 



K 

E 

fc=i 



N 



C{a\q)-Y,Cn{q)Cn{a' 



n=l 



niin (4.20) 



The problem of the inversion is hence reduced to the solution of a set of linear equa- 
tions. It has been shown (Efros et al., 1999; Efros, 1999) that fitting on a grid such 
that K ^ N regularizes the inversion, thus solves the stability problem. 
In order to increase both accuracy and control over the inversion, one can take some 
important measures: 

Physical considerations Often there are physical considerations on the asymptotic 
and threshold of the response function. The choice of a physical shape for /„ 
can diminishes the number of needed expansion functions, while increasing the 
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stability and accuracy of the inversion. Such a set of functions is: 

/„H = c^""e-^. (4.21) 

This choice uses the fact that the threshold behavior is a power law, with 
riQ ~ 1.5, and that the high energy tail decays exponentially with characteristic 
width AE. 

Sum Rules Integration over the entirety of the response with a specific weight func- 
tion depends on ground state properties. Thus, the inversion can be checked 
for consistency. In Chap. 6 an extensive discussion on the subject of sum rules 
in photodisociation processes is given. 

Stability with respect to ai The LIT can be calculated for various widths aj. A 
stability of the inversion with respect to cr/ shows that there are no hidden 
structures of width less than o"/. 



r 



(^ 



Nuclear Wave functions 



The nuclear problem regards the properties of an A-hodj quantum mechanical non- 
relativistic system subjected to strong correlations. As a result, the calculation of 
exact solutions to the problem is a difficult task, demanding elaborate theoretical and 
numerical methods. In order to solve the LIT equations, one needs a good knowledge 
of the ground state wave function, and the ability to calculate matrix elements. 
In the non-relativistic approximation the nuclear Hamiltonian, 

^ = E^+ E^^'^(^"'^^)+ E V^'\r,,r,,n), (5.1) 

i=l i<j=l i<j<k=l 

includes not only a nucleon-nucleon (NN) interaction V^'^\fij), which depends on 
the distance between the nucleons fij = fi — fj, but also three nucleon force (3NF) 
V^^^(fi,rj,rk)- The need in 3NF is observed already in the trinuclei binding energies, 
where state of the art NN potentials lead to underbindings of about 10%. 

In this chapter I will present the use of hyperspherical coordinates and functions 
for the evaluation of nuclear wave functions and matrix elements. The effective in- 
teraction in hyperspherical harmonics (EIHH) approach is introduced at the end of 
the chapter. 

5.1 Few— Body problems in terms of the Hyper- 
spherical Harmonics 

The most simple A-hodj problem is the 2-body problem. One of the conclusions from 
the solution of the 2-body problem is that a natural selection for representing the 
problem is to separate the center of mass movement and the relative motion among 



43 



Chapter 5. Nuclear Wave functions 44 



the bodies. In the case of an isolated system the center of mass movement is trivial 
and does not affect the system. The center of mass coordinate, for a system of A equal 
masses is Rem = -j J2i=i^i- "^^^ internal coordinates are defined as r-=ri — Rem- 
This way, when working in the center of mass frame, Y2^ f- = 0. 

It is clear that a good choice of internal coordinates can help representing the 
nuclear problem. The Jacobi coordinates provide such a choice, as the logic of their 
construction is the separation of a subset of particles from the center of mass of the 
other particles. 

For a one body operator, one chooses 



k-1 



k-l 



l^J^rlj ;A; = 2,...,A (5.2) 

In this set, the last Jacobi coordinate is proportional to the coordinate of the last 






particle, t/a-i = y i^t^a- 

For a two-body operator, one chooses 



A 




Vk = \l^r—. — rA^^-^^Tk >^ ^"O ■,k = l,...,A-l. (5.3) 

In this set, the last Jacobi coordinate is proportional to the relative distance between 
the two last coordinates, tJa-i = \ \ifA-i — "Oi)- 

As a result, the matrix element depends only on the last Jacobi coordinate. One 
can use the invariance of the matrix element to the other Jacobi coordinates by 
defining a generalization of the spherical coordinates, the hyperspherical coordinates. 

5.1.1 Hyperspherical Coordinates 

The hyper-radius of a k-hodj internal system is defined as (Efros, 1972), 

k k+l 

pi=J2^^ = lJ2(^^^--iy- (5-4) 

1=1 i<j 
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From this definition it is clear that p\ = pi_i+Til, so that {pk^i^rjk, Pk) is a pythagorean 
triplet, with a hyper-angle ifk ^ [O, |]: 

Pk-i = pk cos ipk, (5.5) 

r]k = pk sin if k. 

With this, the A — 1 internal degrees of freedom are described by one hyper-radius 
Pa-1, a — 2 hyper-angles (p(A-2) = {^2, ■ ■ ■ , '^A-i}, and 2{A — 1) degrees of freedom 
related to the solid angles of the Jacobi coordinates fi^A-i) = {Vi^ • • • ; Va-i}- 
The volume element related to the internal yl-body system is (Barnea, 1997) 

dV3(A-l) = pT-^i'^PA-idSsA-A (5.6) 

= P^/ji'^pA-1 sin^iipA-l) COS^^~\ipA-l)dipA-ldfjA-ldS'iA-7, 

where dSjq is a hyper-sphere surface element in A^ dimensions. 

One still has to write the kinetic energy operator using hyperspherical coordinates. 
It is again useful to start with an analogy to the 2-body problem, where the Laplace 
operator for one relative coordinate is 

A^, = A,, - 4 (5.7) 

here the radial part is A^^ = -^-^ ( Vi^ ) • Fo^ a set of A — 1 Jacobi coordinates, the 
Laplace operator is a sum A(^_i) = X]i=i A^.. After a transformation to pa-i,^a, 
one gets (Efros, 1972): 



A(A-i) = ^pA-, - -#^ (5.^ 

Pa-i 



with the hyper-radial part 



d^ 3(A-1)-1 d 

dpA-i Pa-1 dpA- 



^PA-^ = TT^ + ^—^ TT^- (5.9) 



The grand-angular-momentum operator Ka-i is built recursively, starting with Ki 
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li. The construction is done by choosing two subsets of Jacobi coordinates, and using 
Eq. (5.5). If one chooses A — 2 coordinates as a subset, the recurrence relation is 



A 



(A-l) 



Ka- 



^{A-2) + ^fjA-i- Putting Eqs. (5.7,5.8) in this relation, one gets: 
92 3p-l)-2] - [3(A-l)-2]cos2(^^_i d 



dip^ 



+ 



A-l 
^A~2 



sin2(y9 



A-l 



dip 



+ 



(5.10) 



A~l 



''A-l 



COS^ (fA-l sm (fA-l 

If we define the internal angular momentum operator La-i = Yli=i ^ii i^ i^ clear 
that the operators Ka-i, Ia-i, Ka-2, -^a-i ^^^ L^_i are commutative. Thus, an 
eigenstate of the Laplacian can be labeled by the quantum numbers Ka-i, ■ ■ ■ , Ki, 



h, 



I 



A-l, 



M-I5 



Li and Ma-i- 



5.1.2 The Hyper spherical Harmonics Functions 

The hyperspherical harmonics functions are the eigen-functions of the Laplace op- 
erator. As explained in the previous section the operator is built recursively, a fact 
which reflects upon the construction of its eigen-functions. 

In general, when treating the nuclear problem as A identical particles (i.e. when 
identifying protons and neutrons in the isospin symmetry) the internal wave func- 
tions belong to an irreducible representation (irrep) of the symmetric group Sa- 
This simplifies the calculation of matrix elements, as for a generic A^-body opera- 



tor O 



{N) 






<i2<...<iN- 



-1 ^n,«2,...,«]v ^^^ problem reduces to a matrix element of an 



operator acting on one A^-body subsystem (^ 



O 



(TV) 



^'> = a)(* 



o 



A-N+1, 



The relevant cases of one-, two-, and three-body operators are written explicit 



y as. 



^ 



0(1) ^') = A{^ Oa ^'> 



{^\0i2)\^') = '^^\ ^\ ^\Oa,a-i\^') 

{^ 0(3) ^'> = ^(A^}M^{^ dA,A-l,A-2 



^') 
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In order to construct a hyperspherical function with such a property, one starts 
in the 2-body problem. In this case there exists one Jacobi coordinate, and the 
hyper-spherical formalism coincides with the usual spherical formalism, thus the 
eigen-f unctions are the spherical harmonics, Yij^miivi)- The latter also belongs to 
a well defined irrep of the symmetry group S2- Adding an additional particle is done 
by the usual coupling of angular momenta, 

^L2,M2;ii,hi^{2)) = Y.mi,m2(^i'^^^2m2\L2,M2)Yi^mi{Vi)Yhm2iV2)- This docs uot ex- 
haust the eigen-function of the 3-body grand-angular momentum K2, as there is an 
additional dependence on the hyper-angle coordinate (p2- The resulting differential 
equation is separable, bringing to the hyperspherical function for three particles, 

3^/f2,-L2,M2,h,/2(^(2),</52) = '^ K2;lihi^2)^L2,M2;li,l2- (^-H) 

By a linear combination of these functions with common A'27 L2, M2 one can construct 
basis functions that belong to well defined irreps of the orthogonal group 0{2). Con- 
sequently, the hyperspherical function depends on the Gel'fand-Zetlin (GZ) symbol 
A2, and a multiplicity label 02- So the hyperspherical function is labeled by five good 
quantum numbers, 3^[a'2](^{2), V^2), with [K2] = i^2, -^2, ^2, A2, 02- 

As mentioned, the construction of the hyper-spherical function for the A-body, 
A — 1 Jacobi coordinates, problem is done recursively. Let us assume the construction 

led to y[KA-2]i^{A~2), V{A~2)) with [Ka-2] = Ka-2, La-2, Ma-2, Ka-3, La-3, Ia-2, Aa-3, aA-3- 

By coupling this to Yi^_,mA^^i'nA-i) one gets 

^L2,M2;h,hi^{2)) = (5.12) 

= X] {LA-2MA-2lA-imA-l\LA~l,MA-l)y[KA-2]i^iA-2),^iA~2))yiA-imA-i{VA-l)- 

MA-2,mA-l 

Next, one constructs the orthonormalized eigen-f unctions of the grand angular mo- 
mentum operator K\_^ (Efros, 1972): 

'^KA-r,lA-lKA-2i.'PA^l) = (5.13) 



1 It- I 3A-8\ 



Ar4-l(AVi; lA^lKA^2)ismifA^iy^-'icOSifA-lf^''PnA-T~''''^'"^ ' \cOs2^ 



A-ly 



Chapter 5. Nuclear Wave functions 48 



where Pn are the Jacobi polynomials, ua-i = ^ ^ — |-^ — ^-^, and the normaliza- 
tion factor is (Efros, 1972): 



Ma~i{Ka~i;Ia~iKa-2) 



{2Ka^i + 3^ - 5)nA-i!r(nA-i + Ka-2 + U-i + ^) 



^ T{nA-i+lA-i + ^)r{nA-i + KA-2 + ^) 

(5.14) 

The function \& is an eigen-f unction oi K\_-^ with an eigen-value -ffA-i(-^A-i+3A— 5). 
By construction, Ka-i — {Ka-2 + ^A-i) > 0, and even. As a result ua-i is a non- 
negative integer. 

The result of the complete separation between hyper-angular coordinate and the 
Jacobi coordinates, means that one can construct the hyper-spherical function 

y[KA-i]i^iA-l), ^{A-l)) = ^KA-i;lA-iKA-2i^A^l)^LA-i,MA-i;[KA-2]lA-ii^{A-2)'f(A-2)) 

(5.15) 
where [Ka-i] stands for Ka-i, La-i, Ma-i, Ka-2, -^a-2, ^a-i, ^a-2, otA-2- These func- 
tions are a complete set, ortho normal in each of the quantum numbers. 

However, using this construction, the hyperspherical harmonics functions still do 
not possess the symmetry property for quantum symmetry, i.e. they do not belong to 
an irrep of the symmetry group Sa- For this purpose we use an algorithm developed 
by Barnea and Novoselsky (1997), which uses an extension of the Novoselsky, Katriel, 
Gilmore (NKG) algorithm (Novoselsky et al., 1988), and yields an A particle state 
belonging to well defined irreps of both the orthogonal and the symmetry groups. 

The process starts with a part overlooked in the previous explanation, the recursive 
construction of the hyperspherical functions as an irrep of 0{A — 1). This step is done 
by understanding that the GZ states can be expressed recursively Atv = [AatAtv-i], 
where X^ is an irrep of 0{N) containing the irreps represented by the GZ state Ajv-i. 
Using this relation, one gets 

JA'ArLTvMTvAAraAr) = (5.16) 

= ^_^ [{KN_iLj^_iX]\r_iaiy_i;lN)KNL]\[\}K]yL]\rX]\[_ia]\[_i]x 

^\{KN^lLN^lXN^lOiN-l]lN)KNLN)KNLNMN). 
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The transformation coefficients connect A^ — 1 coordinates with the A^*^ coordinate, 
thus called the hyperspherical orthogonal coefficients of fractional parentage (hsocfps). 
in Barnea and Novoselsky (1997) it is shown that the hsocfps depend only on the 0{N) 
and 0{N — 1) irreps, and are independent of M^. 

The next step is the reduction of the kinematical symmetry to the permutation 
symmetry, i.e. 0{A — 1) ^ Sa in the case of A identical particles. The NKG 
algorithm is based on the recursive construction of a state of an irrep of Sn from 
a state of an irrep of Sn-i- The algorithm uses the transposition representation 
of a permutation, by which any contained permutation is represented as a subset 
of the transpositions. By construction, a state with A^ particles is characterized 
by the symmetric group irrep that correspond to the group-subgroup chain Si C 
^2 C . . . C Sn, represented by Young diagrams Fi, . . . , T^. For short, one uses the 
Yamanouchi symbol Y^ = [FatFtv-i • • • Fi] = [FatFat^i]. The symmetry adaptation 
of the Oat-i symmetry to include the permutation symmetry has its price, letting 
go of the kinematical symmetry chain of group-subgroup Aat-i and holding only 
the On-1 representation invariant, i.e. Aat-i- Also, a possibility of degeneracy of 
the permutation states in the kinematical group, suggests an additional label (3n- 
This is a kind of parentage transformation, thus its coefficients are called orthogonal 
coefficients of fractional parentage (ocfps). 

Summing this up, we have the recursive relation: 

\KA-iLA-iMA-iXA-iaA-iTAYA-i(3A) = (5.17) 

= Yl [(Aa-2Fa_i/3^_i)A^_i|}A^_iFa/3a]x 

X 2J [{KA-2LA-2>^A-2<yA-2;U-l)KA-lLA-l\}KA-lLA-l><A-l<yA-l]x 

KA-2LA-201A-2IA-1 

X \{KA-2LA-2XA-2aA-2^A-lYA-2(3A-l', 1'A-i)Ka-iLa-iMa-iXa-i) ■ 

Where of course, 
{n\KA-iLA-iMA^iXA^iaA-iTAYA-i(3A) = yKA.,LA-,MA-,XA-.rAYA-.o.'<m (5.18) 
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with Q = {^{A~i),f{A-i)) for the angular coordinates, and a^ = {aA-i, Pa) is the 
multiphcity labeL 

5.1.3 The Hyperspherical Formalism with Internal Degrees 
of Freedom 

Pauh principle suggests that a set of identical particles would have good symmetry 
properties with respect not only to the configurational exchange of two particles, but 
to all the degrees of freedom of the particle. In the case of nucleons, this means that 
the wave function of the A-hodj nucleus can be written as, 

^{Vl ...f]N,Si...SA,ti...tA)='^ RlK]-^{p)H[K]i^, Si . . . sa, ti . . . ^a)- (5.19) 

[K]u 

R[K];u{p) are solutions to the hyper-radial equation, with quantum number z/, and 
[K] = Ja, Ja^ [-^a], Sa, Ta, T4, a^. The function H[k]{^] si . . . SA]ti . . . t^) possesses 
the correct symmetry property for a set of fermions, i.e. antisymmetric in any ex- 
change of particles. The nuclear states are described by total isospin, T, and total 
angular momentum J. The latter is the sum of the total spin S and orbital angular 
momentum L. 

It is clear that in order to couple the hyperspherical harmonics to a spin-isospin 
state, one needs to create the latter as a representation of the permutation group. 

This is done, once again, by means of the NKG algorithm. One uses spin-isospin 
coefficients of fractional parentage (stcfps), to build the state 

XSASXTATXTAYA-iai'^i^i^ ■ ■ ■ , S A]ti, . . . ^Ia) = 2^ (5.20) 

Sa-iTa-io^Y^^ 

[{Sa^i\ sa)Sa{Ta~i\ tA)TATA-ia'^A-i\}SATA^AaT]^ 

[XSA-.TA-.TA-.YA-.aST_^ ® SAtA^^'^-^^^K 

where Sa (Ta) is the total spin (isospin), with projection S\ (T4), and the label a^ 
is needed to remove the degeneracy. 
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Consequently, the function H[k-^{Q; Si, . . . ,SA',ti, . . . , 1^) = {^, Si . . . sa, ti . . Aa 
[K]) is just, 

H[K]in; si, . . . , Sa; h,...,tA) = J2 ^f#T^ (^-21' 



X 



where A-Pa,Ya-i ^^^ phase factors, which in the case of fermionic system are positive 
(negative) when the number of boxes in the diagram Ta below the row of the A*^^ 
particle is odd (even). | F^i | is the dimension of the young diagram. This function 
couples the total spin and orbital angular momentum to angular momentum J a, with 
projection J\, i.e. an LS scheme. However, in microscopic problems, where the forces 
are not necessarily central, a J J coupling scheme is usually preferred. 

In the J J coupling scheme, one first couples the spin and angular momentum of 
the last particle sa = \ and I a to ja, which is then coupled to the total angular 
momentum Ja-i of the residual A — 1 particles, given by the coupling of Sa-i and 
La-2- Finally, the hyperspherical spin-isospin basis state is written as 



AJX 



l[^])(l)=E%r=f ^ [i^A-2TA-l(3U)^A-l\}XA-lTA(3: 

y ^ [{K A-2L A-2XA-2I3 A-2^ Ia-i)Ka-iLa-i\}Ka-iLa-i\a-iPa~i\^ 

Ka-2La-2P^_iI-A-1 

Y, [{SA^i;sA)SA{TA-i;tA)TATA-ia'/_,\}SATATAaf]x (5.22) 

{Sa~1 Sa Sa 

La^2 Ia^i La^i > X 
Ja-1 ja Ja 

\{TA-utA)TATX)x 

I {{Sa-I', Ka-2La-2Xa-2/3a-2^ A-iyA-2l3A-l)JA-l'-, {sa'-,Ia-i)Ja) JaJa) 

This is relevant mainly when the choice of Jacobi coordinates is according to Eq. (5.2), 



Chapter 5. Nuclear Wave functions 52 



in which the last coordinate is proportional to the last particle coordinate. 

A different J J coupling scheme is used when the choice of Jacobi coordinates 
is according to Eq. (5.3), i.e. the last coordinate is proportional to the relative 
coordinate of the last two particles. In this scheme, one couples the last two particles 
total spin s = sa + sa-i and orbital angular momentum I = Ia + Ia~i to a total 
angular momentum j. The latter is then coupled to the total angular momentum of 
the residual A — 2 particle system, which is the coupling of its spin 5*^-2 and orbital 
angular momentum La-z- This leads to the following representation of the basis 
state, 

|[^])(2)= E %#=f%#^ E [i>^A^2TA-l(3U)\A^,\}XA-,TAP\]^ 
YA-.YA-2 VI^aI VUA-1 I ;,^^_^^A_^ 

E \{SA-l\SA)SA{TA-x\tA)TA^A-iaTx\}SATA^AO-f\->^ 

J] \{SA-2\S)SA-x{TA-2\t)TA-l^A-'lO.T2WSA-2TA-2^A-20-f_^-X (5.23) 

SA-2TA-20L^A-1^t 

Y^ (-)— -V(25a-i + 1)(2s + 1)( ^\' \ ^'-' \ X 

jstJA-2 I 2 ^A ^ J 

(-)^-^+^V(2T^-i + l)(2t + l)| ^X' J^ ^^-^' I X 

{Sa~2 s Sa 
La-2 I La-i > X 
Ja-2 3 J A 

\{TA-2,t)TATX)x 

I {{Sa-2',Ka-2La-2Xa-2(3a-2^a-iYa-2(3a-i)Ja-2', {s;l)j) JaJa) 

With these results, one can now proceed to the calculation of matrix elements. 
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5.1.4 Matrix Elements within the Hyperspherical FormaUsm 

In nuclear problems one has to consider matrix elements of operators of pure isospin 
T, and its projection 7^, and spin-orbital coupling to a total angular momentum J 
with projection M.. This operator can act either on a single nucleon, or more. In 
the general case of an operator acting on A^ nucleons, one considers the following 
operator: 

A 

^{N) ~ / ^ ^ V ii) • • • ) ''"ijv) -Sii, • • • 5 Sjjy jC [tij^, . . . ,ti^) [0.Z4:) 

*l<i2<---<«iV = l 

In the case of A^ = 1 or 2, i.e. one and two body matrix elements, one can use hyper- 
spherical wave functions of Eqs. (5.22, 5.23), which lead to explicit dependence on one 
Jacobi coordinate. The advantage is that usually one can separate the hyper-radial 
part from the hyper angular part, i.e. 

QJM ^ n^(Lr]A-i)0^^{7]A_i, spins) (5.25) 



where L = v2 for two-body operator, and L = d ^zi for one-body operator. The 
matrix element of the scalar operator TZ'^ [LrjA-i) includes two integrations. The first 
integration is over the hyper-angle, 

M^{p) = {Ka-21;Ka-i I n-^iLpsinifA) \ KA-2I'; K'a_,) = (5.26) 

Af{KA-i; IKa-2)M{K'a_^-, I'Ka-2) x 

/•f 
X / rf(/;(sinv5)"+°'(cos(^)^~^P^°'^)(cos2y?)7^^(Lpsiny?^)P^"''^)(cos2y?) 
Jo 

where a = / + |, a' = /' + |, and (3 = Ka~2 + ^^y^ (^^^ -^^- (5-13) for the other nota- 
tions). This integration can be performed analytically by expanding the interaction 
to power series (Barnea, 1997). Due to the orthonormality of the Jacobi polynomials, 
this integration includes a finite number of elements. 

One now integrates over the hyperradius, by expanding the hyperradial part of 
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the wave function (cf. Eq. (5.19)) as, 

flw»(p) = Ec£"'^"<'"(9- (5-27) 

Here, 0n"' (x) are the ortho-normahzed generahzed Laguerre polynomials of order Up 
and parameters a, b. The parameter 6 is a length scale for the integration. Explicitly, 
we take 



{R^^yM I M^{p) I rikvAp)) = Ecif "'c-ij'^" r<'w'-*-vir'" (9 M^w':, ^^ 



As a result, the hyperradial matrix element is written as 

(5.29) 
This integration is done numerically, using Gauss-Laguerre quadrature. 

In the current work, we confront also three-body matrix elements, which include 
the coupling of two Jacobi coordinates (Barnea et al., 2004). 

5.2 Effective Interaction with the Hyperspherical 
Harmonics 

The Hyperspherical Harmonics expansion has the correct asymptotic behavior, thus 
provides a promising channel to the solution of the nuclear problem defined by the 
Hamiltonian in Eq. (5.1). Two of the advantages in using this expansion are the 
removal of the center of mass coordinate; and the localization of the quantum de- 
scription, provided by the explicit dependence in the hyperradius, which holds within 
information regarding the entire nucleus. 

Alas, in order to reach convergence in the expansion of the nuclear wave function 
one has to expand up to very high grand angular momentum K. In order to cope 
with this disadvantage it was proposed to replace the bare potential with an effective 
interaction, which truncates the model space. The concept of effective interaction has 
been used traditionally in the framework of harmonic oscillator basis. However, lately 
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(Barnea et al., 2000, 2001b, 2004) it was extended to the hyperspherical harmonics 
basis. This method, the Effective Interaction with the Hyperspherical Harmonics 
(EIHH), has been shown to give excellent results both in comparison to other few- 
body methods, and in confronting previously unreachable nuclear reactions. In this 
section I will present a brief overview of the method, as used for the calculations 
shown in the following chapters. 

5.2.1 Lee— Suzuki Approach for Effective Interactions 

The effective interactions are projections of the bare interactions to a subspace of the 
entire model space. Defining P as the projection operator onto the model space, then 
Q = 1 — P is the projection operator onto the complementary space. By definition, 
the effective interaction leaves unchanged the low lying eigen-energies. 

The first step in the construction of the effective interaction is choosing the model 
space. It is helpful to note that the hyper-radius, by definition, limits the inter- 
nucleon distance. This fact makes it possible to choose the model space as the product 
of the entire hyper-radial space and the set of HH basis functions with K < K^ax- 
It is clear that V^jf — > V^ as K^ax -^ oo, and the same follows for the eigen-energies 
and states of the Hamiltonian. 

In order to continue, one method is to use the Lee-Suzuki procedure (Suzuki and 
Lee, 1980, 1982, 1983), in which one uses a similarity transformation X. The effective 
Hamiltonian is given by 

Ueff = PX-^nXP. (5.30) 

The effective interaction acts only in the P-space. The operator X = e'^ re-projects 
the information contained in the Q-space back to the model space, i.e. u = QujP 
(which also leads to u;^ = u;^ = . . . = 0, and thus X = 1 +uj, and X~^ = 1 — uj). One 
can summarize this method in the following way (we use the notation Oab = AOB) 
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(Bacca, 2004): 
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(5.31) 
The implicit condition that u has to fulfill Ti-Qp = can be rearranged to the equation: 



Q{n + [n,ij] - ijnuj)p = 0. 



(5.32) 



The strength of the effective interaction approach is that the low energy solutions 
to the original problem can be exhausted from the eigen-vectors of the effective 
Hamiltonian, 'Heffi which by construction are in the P space. The theorem states: 
let Ip be the eigen-vectors of this subspace, then ep = Xip are eigen-vectors of the 
true Hamiltonian with the same eigen-value. 

In order to prove this theorem, let | \&^) G ep, i.e, Heff \ ^^^) = e^ \ ^^) and 
P I ^^> =1 ^^). Then n\^^)=nX \ ^^) = XX-^nX I ^^) = XHeff I ^^> = 
e^X I \&^) = e^ I "^fj.), which concludes the proof. 

The subspace ep contains information regarding both the P and the Q spaces. 
Furthermore, the similarity transformation depends only on this space. Let | a) 
(I /3)) be a basis for the P (Q) space, and {A)a^ = {a \ \I/^), {B)^f^ = (/3 | \i/^). 
Then, it is easily proven that a solution to Eq. (5.32) is 



UJ 



A. {BY 



(5.33) 



This similarity transformation is not hermitian, however an equivalent effective 
interaction can be written as, 
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eff 
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P{1+UJ^U)P ^P{1+UJ^U)P 



(5.34) 



It is clear that the effective interaction stores information regarding correlations 
in the Q-space, while solving only in the P space. Alas, this procedure produces, in 
general, A-hodj interactions, which solution is as difficult as finding the full space 
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solutions. To cope with this problem, one approximates the effective interaction, 
while keeping the property of converging to the bare interaction as P — *> 1. 

A good approximation is used in the "No Core Shell Mode" (NCSM) approach 
(Navratil and Barrett, 1996, 1998, 1999). The approximation is to include in the 
effective potential only two-body effective interactions: 



V.St = PX'' 






xp -^ v:ff = y: y%r (5.35) 



The two-body effective potential is the non-trivial part of the two-body effective 
Hamiltonian, i.e. one decomposes HS^'^ = ^o + ^i^j=i Vij with [Hq, P] = [Hq, Q] = 
and QHqP = PHqQ = 0. Then the two-body effective Hamiltonian is built using a 
Lee-Suzuki transformation: 

nf^f = ^^ + ^^ 7^(^) ^'^""^ =, (5.36) 

thus, 

A 

E ^i/1.- = ^fA - ^HoP. (5.37) 

5.2.2 EIHH method 

In the EIHH method, one solves the hyper-radial part of the Hamiltonian. Then, for 
each hyper-radius p, the two-body hamiltonian is clearly 

Let us note that: 

• This Hamiltonian includes the total hyper-angular kinetic energy. 

• The only interacting particles are the "last" two. 

• The residual A — 2 system acts as a mean-field force, entering only through the 



Chapter 5. Nuclear Wave functions 58 



collective hyper-radius. The hyper-radius enters not only in the kinetic energy, 
but also in the potential matrix element, as the relative distance between the 
last two particles is f^^A-i = V^Va-i = {V^psiT^^A,VA-i)- 

• The calculation includes only one-dimensional integration. 

Choosing the model space within the HH formalism is quite natural, as the two-body 
Hamiltonian of Eq. (5.38) is diagonal in the quantum numbers of the A — 1 residual 
system. Thus, the P2^space is chosen as P2 = {[-f^] : K < K^ax, [Ka-2] fixed}. It is 
clear that Q2 = {[K] : K > K^ax, [Ka-2] fixed}. A numerical calculation of the Q2 
space includes a numerical cutoff Kmax which assures convergence. This depends 
on the orbitals of the potential, and was checked to be Kmax ~ 60 for S-wave 
interaction and Kmax ~ 180 for P-wave interaction. The construction of ep space 
is now trivial, using Eq. (5.33). 

By construction, the effective interaction is "state dependent". This means that 
for every fixed value of Ka-2 one has to calculate a similarity transformation. In 
Mintkevich and Barnea (2004) it was demonstrated that the "no-core" approach is 
equivalent to the bare Hamiltonian (through a unitary transformation), in the limit 
P — > 1. It was also shown that any effective two body operators constructed through 
the Lee-Suzuki similarity transformation of Eq. (5.36) can be regarded accurate to 
second order in the limit P — > 1. 

The EIHH method accelerates the convergence of the HH expansion, and has been 
shown to do so in the calculation of nuclear bound states and reactions of nuclei in 
the mass range 3 < A < 7 (for example (Barnea et al., 2000, 2001b, a; Bacca et al., 
2004, 2002; Gazit and Barnea, 2004; Gazit et al, 2006a; Gazit and Barnea, 2007c)). 




Photoabsorption on ^He 

li 

WITH A REALISTIC NUCLEAR 



FORCE 



In previous chapters I outlined a detailed description of the theoretical background for 
ab-initio calculation of few-body nuclear structure and reactions. In the first section 
of this chapter, the application of these methods to the ^He total photoabsorption 
cross section is presented with the realistic nucleon-nucleon (NN) potential Argonne 
fi8 and the three- nucleon force (3NF) Urbana IX, as published in Gazit et al. (2006a). 
In the second part of the chapter, based on (Gazit et al., 2006b), sum-rules of this 
cross-section are investigated, emphasizing relations to ground state properties. 

6.1 Photoabsorption Calculation 

6.1.1 The Unretarded Dipole Approximation 

For low energy reactions, below the pion production threshold, a good approximation 
for the photoabsorption cross-section, a^, is using only the leading electric dipole 
{El) in its non-relativistic long wavelength approximation. In addition, since the 
Siegert theorem (Sec. 2.3) applies for the conserved electro-magnetic current, the 
dipole operator is proportional to the Coulomb multipole. Thus, the effects of MEG 
are implicitly included. 

As a result, the total photoabsorption cross section is given by 

^tH = 2J + ^ ^^i^) ' (6-1) 

where a is the fine structure constant, Jq denoting the nucleus total angular momen- 
tum, and 

Ri.)^f:\i,AD.\%)fm-E.-.) (6.2) 



59 



Chapter 6. Photoabsorption on "^He with a realistic nuclear force 60 



^3,/ 



is the response function in the unretarded dipole approximation with D^ = X]i=i "V^- 
The wave functions of the ground and final states are denoted by |^&o//) and the 
energies by -Eo//; respectively. The operators rf and z^ are the third components of 
the i-th nucleon isospin and center of mass frame position. 

We refer back to Sec. 2.3 and recall the Siegert theorem, which ensures that the 
dominant part of the exchange current contribution is included. In the classical few- 
body nuclei this has proven to be an excellent approximation (Golak et al., 2002; 
Arenhovel and Sanzone, 1991), particularly for photon energies below 50 MeV. 
For triton, cr^(^H) (Golak et al., 2002) the contributions of retardation and other 
multipoles lead to an enhancement of a^ by less than 1% for u < 40 MeV. The 
contribution grows with energy transfer. This result is rather insensitive to the nuclear 
force model. When calculating with AV18 NN potential only, the contribution is 5% 
for 00 = 60 MeV, 16% for uj = 100 MeV, and 26% at pion threshold {uj = 140 MeV). 
When including UIX 3NF the numbers change to 5%, 18%, and 33% at uj = 60, 
100, and 140 MeV, respectively. If one bares in mind that the cross-section is almost 
exhausted in the giant dipole resonance (GDR), below 40 MeV, the approximation is 
appropriate for the photoabsorption process below pion-threshold. 

6.1.2 Calculation of The Response Function 

In order to find the Response function, the LIT method is used, calculated with the 
Lanczos algorithm (see Chap. 4). 

The EIHH expansions of the ground state |\l/o) and the LIT vector |\1/) are per- 
formed with the full HH set up to maximal values of the HH grand-angular momentum 
quantum number K {K < K^ for ^q, K < K^ for \1/). 

The ground state properties of ^He calculated using the EIHH approach are shown 
in Tab. 6.1. The convergence of the binding energy, Eb, and matter radius are pre- 
sented as a function of K^^. Since the EIHH method is not variational the asymptotic 
Eh value can be reached from below or above, in fact both cases are realized in Tab. 6.1. 
Both convergence patterns are not sufficiently regular to allow safe extrapolations to 
asymptotic values. The calculated ground state properties of ^He agree quite well 
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Table 6.1: Convergence of HH expansion for the ^He binding energy E^ [MeV] and 
root mean square matter radius (r^)2 [fm] with the AV18 and AV18+UIX potentials. 
Also presented are results of other methods (see text). 





AVI 8 




AV18+UIX 




K^ 

-'^m 


E, 


(r^)^ 


Eb 


{rr^ 


6 


25.312 


1.506 


26.23 


1.456 


8 


25.000 


1.509 


27.63 


1.428 


10 


24.443 


1.520 


27.861 


1.428 


12 


24.492 


1.518 


28.261 


1.427 


14 


24.350 


1.518 


28.324 


1.428 


16 


24.315 


1.518 


28.397 


1.430 


18 


24.273 


1.518 


28.396 


1.431 


20 


24.268 


1.518 


28.418 


1.432 


FYl 


24.25 


- 


28.50 


- 


FY2 


24.22 


1.516 


- 


- 


HH 


24.21 


1.512 


28.46 


1.428 


CFMC 


- 


- 


28.34 


1.44 



with those of other methods, presented in Tab. 6.1. The other calculations are done 
with the FY equations (in Nogga et al. (2002); Lazauskas and Carbonell (2004)), 
the HH expansion (in Viviani et al. (2005)) and the GFMC method (Wiringa et al., 
2000), and appear in this order in Tab. 6.1. One should compare the binding energies 
to the experimental value of 28.296 MeV. 

The EIHH convergence of the transform C{a) is excellent for the AV18 potential. 
Here we discuss in detail only the case AV18+UIX, where the convergence is quite 
good, but not at such an excellent level. The reason is that in our present EIHH 
calculation an effective interaction is constructed only for the NN potential, while the 
3NF is taken into account as bare interaction. 

In Fig. 6.1 we show results for the transform C obtained with various Km and K^ 
values. Since we take K^ = Km — 1, the corresponding transform can be denoted by 
Cxm- One sees that there is a very good convergence beyond the peak, that the peak 
height is very well established, but that the peak position is not yet fully converged 



Chapter 6. Photoabsorption on "^He with a realistic nuclear force 62 



0.005 



0.004 



Fig. 1 



> 

s 



I 



0.003 - 



2 0.002 

II 



0.001 - 



0.000 




40 
Or [MeV] 



Figure 6.1: Convergence of Cxm with a/ = 10 MeV (AV18+UIX). 



in view of the fact that we aim for a percentage level accuracy. In fact with increasing 
Km the peak is slightly shifted towards lower ctr. 

By observing the relative error of the transform, A^^ = Cxm^w/'^iQ "with Ca,i3 = 
Ca — Cp, this point is illustrated better, see Fig. 6.2 (the chosen a/j-range starts at the 
^He(7) break-up threshold). One again notes the very good convergence for cjr > 30 
MeV with almost identical results from £13 to £19. Altogether we consider our result 
for cr^ > 30 MeV as completely sufficient. However, it is obvious that convergence is 
not entirely reached for lower a^. 

A trivial solution for this is an additional increase in Km- However, £19 includes 
already 364000 states in the HH expansion. Thus, a further increase is beyond our 
present computational capabilities. Fortunately a closer inspection of Fig. 6.2 shows 
that the convergence proceeds with a rather regular pattern: (i) £13,11 ~ £11,9 and 
£17,15 ^ £15,13 and (ii) £19,17 — £17,15/1.5 ~ £13,11/(1. 5)^. Therefore it is possible to 
obtain an extrapolated asymptotic result. 
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Figure 6.2: Convergence of /^Km = {^k^ - l^\<i)ll^\^ (AV18+UIX). 

The extrapolation used is a Pade approximation (Fabian and Arenhovel, 1976; 

Basdevant, 1972) 

^x„ = ^/^.-B + /:x„-4,i.„-8/(l - f^""^"'-' ) • (6.3) 

6.1.3 Results and Discussion 

In Fig. 6.3 the resuhs for u^^k^, obtained from the inversions of the transforms Ck^ 
are presented. Due to the Lorentz kernel the a^ presents the same features as £ itself: 
stable peak height with a value very close to 3 mb, and not yet completely converged 
peak position. 

In Fig. 6.3 the asymptotic Pade approximations are given, where cr^^^^ and cr^^^g 
result from the inversions of C^ and £^. Unquestionably, the extrapolated C^^ have 
a lower numerical quality than the calculated Ck^ and consequently the stability 
for the corresponding inversions worsens. Therefore, an additional constraint in the 
inversion is imposed by fixing the peak cross section to the already converged value of 
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3 mb. In Fig. 6.3 it is evident that a^u and cr^ig are very similar, hence establishing 
a very good approximation for the asymptotic a^. One also notices that compared to 
cr^^ig they show a shift of the peak position by about 1 MeV towards lower energy. 

The final results for the cross-section are presented in Fig. 6.4a. Due to the 3NF 
one observes a reduction of the peak height by about 6% and a shift of the peak 
position by about 1 MeV towards higher energy. Very large effects of the 3NF are 
found above 50 MeV with an enhancement of a^ by e.g. 18, 25, and 35% aX u = 60, 
100, and 140 MeV, respectively. The discussion in the previous section indicates that 
the 3NF effect could change somewhat if all multipole contributions are considered. 

It is very interesting to compare the 3NF effects on o"^(^He) to those found for 
cr^(^H/^He) (Efros et al., 2000; Golak et al., 2002). Surprisingly, the reduction of the 
peak height is smaller for ^He. For ^H/^He the size of the reduction is similar to the 
increase of Ei, (10%), whereas for ^He the 3NF increases Eh by 17%, but reduces the 
peak by only 6% and thus cannot be interpreted as a simple binding effect. 
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Figure 6.4: Total '^He photoabsorption cross section: (a) o^ (MT,AV18) and a^^^g 
(AV18+UIX); experimental data from (Arkatov et al., 1979). (b) as (a) but also 
included upper/lower bounds and various experimental data (see text), area between 
dotted lines (Berman et al., 1980; Feldman et al., 1990), dotted box (Wells et al., 
1992), squares (Nilsson et al., 2005), and circles (Sliima et al., 2005). 
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One can find additional differences at higher energy-transfer. The enhancement of 
(T^('^He) due to the 3NF is significantly stronger, namely about two times larger than 
for the three-body case. Interestingly this reflects the above mentioned different ratios 
between triplets and pairs in three- and four-body systems. An integration of a^ up 
to pion threshold yields 93.0 (AV18) and 97.4 MeVmb (AV18+UIX) thus showing a 
rather small difference. If one considers also contributions beyond pion threshold the 
difference becomes more pronounced, in fact integrating up to 300 MeV one finds 103 
(AV18) and 114 MeVmb (AV18-I-UIX). The latter values correspond to enhancements 
of 72 % (AVI 8) and 91 % of the classical Thomas- Reiche-Kuhn sum rule. Different 
sum-rules will be discussed in the next section. 

In Fig. 6.4a we compare the cross-section with the one achieved using the semi- 
realistic Malfiiet-Tjon (MT) NN potential (Efros et al., 1997b; Barnea et al, 2001a). 
Similar to a^('^H/'^He) (Efros et al., 2000) one finds a rather realistic result in the 
giant resonance region (overestimation of the peak by about 10-15%) and quite a 
correct result for the the peak position; however, at higher energy a^ is strongly 
underestimated for ^He, by a factor of three at pion threshold. Also shown are data 
from Arkatov et al. (1979), the only measurement of cr^('^He) in the whole energy 
range up to pion threshold. In the peak region the data agree best with the MT 
potential, while for the high-energy tail one finds the best agreement with the AV18 
potential. 

A comparison of the low-energy results to further data is presented in Fig. 6.4b. 
For the AV18+UIX case upper/lower bounds are included to account for possible 
errors in the extrapolation, Eq. (6.3). As bounds we take ±(cr?°]^g — (T^^i9)/2; we 
believe that this is a rather safe error estimate. 

For a better understanding of the data in Fig. 6.4b a few comments are needed: 

1. In Wells et al. (1992) the peak cross section is determined from Compton scat- 
tering via dispersion relations. 

2. The dashed curve corresponds to the sum of cross sections for (7,n) from 
Berman et al. (1980) and (7,p)^H from Feldman et al. (1990) as already shown 
in Efros et al. (1997b). 
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3. The data from the above mentioned recent (7,n) experiment (Nilsson et ah, 
2005) are included only up to about the three-body break-up threshold, where 
one can rather safely assume that a^ ^ 2(T(7,n) (see also Quaglioni et al. 
(2004)). 

4. In Shima et al. (2005) all open channels are considered. One sees that the 
various experimental a^ are quite different, exhibiting maximal deviations of 
about a factor of two. 

The theoretical a^ agrees quite well with the low-energy data of (Berman et al., 1980; 
Feldman et al., 1990). In the peak region, however, the situation is very unclear. 
There is a rather good agreement between the theoretical a^ and the combined data 
of (Nilsson et al., 2005) and (Wells et al., 1992), while those of (Berman et al., 1980; 
Feldman et al., 1990) are quite a bit lower. Very large discrepancies are found in com- 
parison to the recent data of Shima et al. (2005). It is evident that the experimental 
situation is rather unsatisfactory and further improvement is urgently needed. 

6.2 Photonuclear Sum Rules and the Tetrahedral 
Configuration of ^He 

Sum rules (SR) are related to moments of different order of the photonuclear cross 
section and reflect important electromagnetic properties of nuclei. In fact they can 
often be expressed in terms of simple ground state properties in a model independent 
or quasi model independent way. Well known examples are the Thomas- Reiche-Kuhn 
(TRK) sum rule (Landerburg and Reiche, 1923; Kuhn, 1923; Thomas, 1925), which 
gives information about the importance of exchange effects in nuclear dynamics via 
the so-called TRK enhancement factor k^^^, the bremsstrahlung sum rule (BSR) 
(Levinger and Bethe, 1950; Brink, 1957; Foldy, 1957; Dellafiore and Brink, 1977), 
which is connected to the nuclear charge radius and to the mean distance between 
protons (Dellafiore and Lipparini, 1982), and the polarizability sum rule (PSR) (Friar, 
1975), related to the electric nuclear polarizability, which is proportional to the shift 
in energy levels due to external electric field. 
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These SR assume that the dominant contribution to the cross section comes 
from unretarded electric dipole (ElUR) transitions. Two- and three-body stud- 
ies (Arenhovel and Sanzone, 1991; Golak et ah, 2002) indeed confirm that other 
contributions are much smaller. Much discussed is also the Gerasimov-Drell-Hearn 
(GDH) sum rule (Drell and Hearn, 1966; Gerasimov, 1965), which is related to the 
nuclear anomalous magnetic moment. 

In this chapter, the TRK, BSR and PSR of '^He are considered within a real- 
istic nuclear potential model consisting of two- and three-body forces (AV18 and 
UIX (Wiringa et al., 1995; Pudliner et al., 1997)). The GDH sum rule is trivial for 
^He: it vanishes, since the '^He total angular momentum is equal to zero. Also inves- 
tigated are the related moments by integrating explicitly the properly weighted total 
photoabsorption cross section calculated in the previous section. 

The aim of this study is to show that in some cases sum rules can allow an exper- 
imental access to two-body properties of the nuclear ground state, like the proton- 
proton, neutron-neutron and proton-neutron distances. In the case of ^He this allows 
to test the validity of the configuration tetrahedral symmetry of this nucleus and at 
the same time to "measure" the amount of symmetry breaking. 

This section aims also at providing a guideline for experiments, for which a direct 
test of SR is difficult, as only lower bounds for the SR can be determined, as well as 
at giving an idea of the reliability of the SR approach to heavier systems, where the 
direct theoretical determination of the cross section, and therefore its integration, is 
presently out of reach. The advantage to perform this kind of study in ^He, com- 
pared to analogous ones in the two- (Arenhovel and Fabian, 1977) and three-body 
systems (Efros et al., 1997a), is that ^He is a rather dense and compact nucleus, re- 
sembling heavier systems more closely. Only now that realistic theoretical results for 
the photonuclear cross section are available such a study is possible and one can put 
the extrapolation of the results to heavier systems on safer grounds. 
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6.2.1 SR and the Nucleus Ground State 

We start by recalling the formalism of the photonuclear SR. The various moments of 
the photonuclear cross section are defined as 

mn{oj) = r rf^o^^afu^H, (6.4) 

where u is the photon energy and Uth and uj indicate threshold energy and upper 
integration limit, respectively. Where a^^^^ioj) is the unretarded dipole cross section, 
which can be written as 



afu^M = ^c.'i?M, (6.5) 



7 

where Q = 47r^a/(3(2 Jq + 1)). R{u!) is the nuclear response function to excitations of 
the operator D is the unretarded dipole operator: 

A 

D = J2f,Tf/2, (6.6) 

where A is the number of nucleons and rf and r^ are the third component of the 
isospin operator and the coordinate of the ith particle in the center of mass frame, 
respectively. For n = the SR is the TRK SR, n = — 1 is the BSR, and n = —2 the 
PSR. 

Assuming that a^^^^{uj) converges to zero faster than u;"""^ and applying the 
closure property of the eigenstates of the hamiltonian H, it is straight forward to get 
the following identities for SR with n = 0,-1, —2: 

S™^ = mo(oo) = |(0| [3, [h,dJ\ |0) (6.7) 

E^^^ = m^i{oo) = g {0\D ■ D\0) (6.8) 

E^SR ^ m.2{oo) = gJ2iEn-Eo)-'\{n\D\0)\\ (6.9) 

n 

Working out the expressions in Eqs. (6.7 - 6.9), one finds that those moments 
are related to interesting properties of the system under consideration. In fact the 
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TRK sum rule is also given by the well known relation (Landerburg and Reiclie, 1923; 
Kuhn, 1923; Thomas, 1925) 

E™^ = ^^ (1 + z.™) , (6.10) 

2mA ' 

where N and Z are the neutron and proton numbers, respectively, m is the nucleon 
mass and k^^^ is the so-called TRK enhancement factor defined as 

«:™-^(0|[A[^,^]]|0). (6.11) 

From this expression it is evident that h^^^ embodies the exchange effects of the 
nuclear potential V ^ since the double commutator in (6.11) vanishes for systems like 
atoms, where no exchange effects are present. 

In the literature one finds a few interesting equivalences for the bremsstrahlung 
sum rule. Rewriting the dipole operator as D = {NZ/A)Rpmi where RpN denotes 
the distance between the proton and neutron centers of mass, one has (Brink, 1957) 

S^'^=^(^)'(0|i?|^|0). (6.12) 



Foldy (1957) demonstrated that 



NZ 
A-l 



J.BSR ^^__^^2^^ ^g_^3^ 



where (r^) is the mean square (m.s.) point proton radius 



i=l 

However, this relation is valid only under the assumption that the ground state wave 
function is symmetric in the space coordinates of the nucleons. 

Dellafiore and Brink (1977) have found that, in the framework of the oscillator 
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shell model, one has 

E^^^=g[z^{rl)-Z{rl^)), (6.15) 

where {r^) is the m.s. distance of protons with respect to the proton center of mass 
Rp 

iO^^mf^i^^-Rpm. (6.16) 

Later, in Dellafiore and Lipparini (1982), it was shown that the validity of Eq. (6.15) 
is not limited to the oscillator shell model, but it is a model independent relation, 
which can also be written as 



S^"^=^^^^(rJ)-^^^^^;r^(r^p)), (6.17) 

where (r^ ) is the m.s. proton-proton distance 



-^--(0|X;(r1-r,)^|0). (6.18) 



For the BSR two additional relations exist, which are easy to prove, but, to our 
knowledge, have not been considered in the literature, viz 



.BSR n I ,r2;2, N{N-1)^, 



S"^" = Q [N^iri) - "^ ' {rD) (6.19) 

and 

S^^^=^^((rJ„)-(rJ)-(0), (6.20) 

where (r^) is the m.s. point neutron radius and (r'^a) are the m.s. nucleon-nucleon 
(NN) distances, i.e. 



1 ^ 

i^L) = ^(^_i) (OlE(^"'--^"':'-)^|Q)' (6.21) 

-. Z N 

(^l) = T77(0|EE(^"'^-^"'^-)'|0)- (6.22) 



NZ 

i=i j=i 
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It is interesting to note that Eqs. (6.17), (6.19) and (6.20) express S^^^ via a one-body 
((r^)) as well as a two-body quantity ((r^/j)). 

Finally, regarding the polar izability sum rule, one has 

S^^^ = 27r2a,,, (6.23) 

where an denotes the nuclear polarizability in the ElUR approximation. 

6.2.2 Calculation of SR using the LIT method 

There are various approaches for the calculation of moments and sum rules. The 
obvious one is to obtain the moments by integrating the ^He total photoabsorption 
cross section of Sec. 6.1. However, a more direct approach to obtain the SR exist, as 
explained in the following. 

One way to evaluate the LIT is by using the Lanczos technique described in 
Subsec. 4.2.1. We recall that the LIT can be re-expressed as 

C{aR, ai) = — {Q\D ■ D|0) Imxoo(^) (6.24) 

with z = Eq + aR + icTi and 

— * 

100) = ^==^. (6.25) 

^J{0\D-D\0) 

Xqq{z) can be expressed (cf. Eq. (4.16) as a continued fraction containing the Lanczos 
coefficients a, and 6,, defined in Eq. (4.10). 

Therefore, the implementation of the Lanczos algorithm leads to C{aR, aj). While 
the inversion of the LIT (Efros et al., 1999) gives access to R{uj), and thus to the 
moments of Eq. (6.4), the normalization of the Lanczos "pivot" \LP) = D\0) and the 



Chapter 6. Photoabsorption on "^He with a realistic nuclear force 73 



Lanczos coefficients allow to obtain the SR of Eqs. (6.7-6.9). In fact one has: 

S^'"' = GxooiEo), (6.26) 

S^^SR = g{LP\LP), (6.27) 

S™K = (ao-^o)S^'''. (6.28) 

6.2.3 Results and Discussion 

We use a HH basis, therefore the ground state |0), the Lanczos "pivot" \LP) and the 
Lanczos coefficients a^ are given in terms of HH expansions. While for the ground 
state the expansion is characterized by an even hyperspherical grand-angular quantum 
number K and total isospin T=0, \LP) has to be expanded on {K' = K + 1,T = 1} 
states (we neglect the AV18 isospin mixing, which is very small, as shown for the ^He 
ground state in Nogga et al. (2002); Viviani et al. (2005)). The rate of convergence 
of the various SR results from Eqs. (6.26-6.28) is given in Tab. 6.2 as a function of the 
hyperspherical grand-angular quantum number K. One observes sufficiently good 
convergence patterns for S"""^^ and S^^^. For the latter an additional test of the 
convergence is performed, calculating S^^^ directly as mean value of the operator 
D ■ D on the ground state (Eq. (6.8)). In this way the expansion of \LP) on {K' = 
K + 1,T = 1} states is avoided. We obtain practically identical results. 

From Tab. 6.2 one sees that the convergence of S^^^ is slower when the three- 
nucleon force (3NF) is included. This is not a surprise, in view of the problem found 
for the cross section itself in the previous section. In fact it has been shown that 
the peak of the giant dipole resonance is slightly shifting towards lower energies with 
increasing K. The sum rule S^^^, which has the strongest inverse energy weighting, 
is more sensitive to this shift than the other two SR. 

We have used cr°° extrapolated in the previous section, to determine from Eq. (6.4) 
the various moments for ui = 300 and 135 MeV. These results are also listed in 
Table 6.2. 

For S"""^^ one sees that the SR is not yet exhausted at 300 MeV. In fact more 
than 20% of the strength is still missing. At pion threshold, one has only about 2/3 
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AV18+UIX 




5.PSU 


J^BSK 


5.TKK 


K 


10-2 [mb MeV-i 


mb 


102 [mb MeV] 


8 


6.230 


2.398 


1.430 


10 


6.277 


2.396 


1.448 


12 


6.331 


2.394 


1.451 


14 


6.382 


2.401 


1.458 


16 


6.434 


2.406 


1.460 


18 


6.473 


2.410 


1.462 




AV18 




7.681 


2.696 


1.383 




AV18+UIX 


u) [MeV] 


m_2 (d)) 


m_i {u) 


mo(cj) 


135 


6.55 


2.27 


.944 


300 


6.55 


2.37 


1.14 



Table 6.2: Convergence in K of the SR for AV18+UIX potentials. The converged 
AV18 results are also shown. The last two lines of the table show the convergence in 
uj of the various moments. 



of S"""^^. As was already discussed in (Efros et al., 2000) for the triton case, the 
rather strong contribution from higher energies seems to be connected to the strong 
short range repulsion of the AV18 potential. As to the TRK enhancement factor 
the present calculation gives k™^ = 1.31 for AV18 and 1.44 for AV18+UIX. These 
numbers are somewhat larger than older results obtained either with a variational 
wave function and AV14+UVII potential (k™^ = 1.29 (Schiavilla et al., 1987)) or 
with more approximated wave functions and various soft and hard core NN potentials 



K 



TRK 



0.9 - 1.30 (Weng et al, 1973; Heinze et al., 1978; Gari et al., 1978)). 



For S^^^, as expected, the contribution at high energy is much smaller than for 
T?"^^ , in fact at cj = 300 MeV the missing sum rule strength is less than 2%. For 
^PSR ^YiQ strength beyond 300 MeV is even more negligible. Actually in this case 
the explicit integration leads to an even higher result than the sum rule evaluation 
of Eq. (6.27). The seeming contradiction is explained by the already discussed fact 
that TF^^ is not yet convergent, while for the explicit integration an extrapolated 
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cross section is used. Indeed, integrating the K=18 cross section we obtain a value 
of 6.46 mb, which is consistent with the corresponding sum rule result of 6.47 mb, 
within the numerical error of the calculation. For the AV18+UIX force the value 
of the polarizability ao that we deduce from the extrapolated S^^^ is 0.0655 fm'^. 
The AV18 result, which already shows a good convergence for K = 16, is 0.0768 
fm^. This means that the 3NF reduces the polarizability by 15%. It would be very 
interesting to measure this nuclear polarizability by Compton scattering, as a test of 
the importance of the three-body force on such a classical low-energy observable. We 
find that S^^^ is the SR that is affected most by the 3NF. In fact E^^^ is reduced by 
only 10% and one has an opposite effect on S^^^ with a 5% increase. The quenching 
or enhancement of SR due to the 3NF is the reflection of its effects on the cross section 
i.e a decrease of the peak and an increase of the tail. 

This work allows a few conclusions regarding a very old question, already discussed 
in Levinger and Bethe (1950), of the "existence" of the SR (finiteness of mn{uj) for 
uj -^ oo), which is connected to the high-energy fall-off of the ElUR cross section. 
Since there is a rather good consistency between the SR and the moment values, it 
can be stated with a rather high degree of confidence that S"""^^ (and consequently 
the other two SR) "exist" . Therefore one can try to extract some information about 
the high energy behavior of the cross section and hence about the "existence" of 
SR with higher n. With an uj~p ansatz for the fall-off of the cross section above 
pion threshold, and requiring that S"'"^^ is 146 mb MeV (see Table 6.2), one gets a 
rather weak energy fall-off, i.e. p ~ 1.5. This value is also consistent with S^^^. In 
fact, adding such a tail contribution to m_i(135) one gets E^^^ = 2.39 mb, to be 
compared with 2.41 mb in Table 6.2. The value of p might be somewhat different for 
other potentials, but probably it will not change much. An additional, rather safe, 
conclusion is that higher order SR do not exist for realistic nuclear potential models. 

A very interesting aspect can be viewed by observing the BSR. As already men- 
tioned S^^^ contains information about one- and two-body densities via (r^) and 
(r'^a), respectively. This means that a measurement of S^^^ and the knowledge of the 
experimental m.s. radius allow to determine (r^^) via Eqs. (6. 17), (6. 19) and (6.20). 
In this way one gets information about the internal configuration of ^He as it is 
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explained in the following. 

In his derivation of Eq. (6.13), Foldy assumed a totally symmetric ^He spatial wave 
function, which corresponds to a configuration where the four nucleons are located at 
the four vertexes of a tetrahedron. For such a configuration one has (r^) = (r,^) = (r^) 
and (r^p) = (r^„) = (r.^^) with Qt = {rl^p) / {r'^) = 8/3. Foldy's assumption is a 
very good approximation for ^He, but other spatial symmetries (mixed symmetry, 
antisymmetric) are also possible. 

What can he learned from S^^^ with respect to this question? 

For ^He, which is a T=0 system one can safely assume that {r'i) = (r^) = (r^) (isospin 
mixing is tiny (Viviani et al., 2005; Nogga et al., 2002)). Using in Eqs. (6.17), (6.19) 
and (6.20) the AV18+UIX value of (r^) = 2.04 fm^ from Sec. 6.1 (which coincides with 
the experimental one (Borie and Rinker, 1978), corrected for the proton charge radius) 
and S^^^ from Table 6.2, one obtains (r^^) = (r^„) = 5.67 fm^ and (r^^) = 5.34 fm^, 
i.e. two values which differ by about 6%. 

The ratios Qpp(nn) = ('"pp(„,n))/('^^) ^^^ Qnp = {rnp)/{r'^) are not much different 
from the correspondent value Qt of a classical tetrahedral configuration. We obtain 
Qpp = 2.78 and Qnp = 2.62 instead of Qt = 2.67. One notices that Qpp(nn) ~ Qt — 
'^{Qnp — Qt)- This reflects the different numbers of proton-proton and neutron- 
neutron pairs (2) with respect to proton-neutron pairs (4). Using Eq. (6.12) one can 
also derive the distance between the proton and neutron centers of mass. One has 
Rpn = 1.58 fm instead of 1.65 fm for the tetrahedral configuration. 

Notice that with (r^) = 2.04 fm^ the "tetrahedral" BSR would be 2.62 mb, the 
same value that one obtains using Eq. (6.13). This value is 9% larger than our result. 
The distortions that I find from the 9% smaller BSR are the consequence of the 
different effects of the potential on isospin triplet and isospin singlet pairs. 

A rather intriguing conclusion arises: when considered in its body frame ^He 
should look like a slightly deformed tetrahedron. Of course this statement has to 
be interpreted in a quantum mechanical sense, regarding the mean square values of 
the nucleon-nucleon distances on the two-body density. It is clear that one cannot 
measure this deformation " directly" , since it is not a deformation of the one-body 
charge density (the ^He charge density has only a monopole). On the other hand such 
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a deformation is accessible experimentally in an indirect way via the measurements 
of the charge radius and of S^^^. 

This leads to the question how exactly S^^^ can be measured in a photonuclear 
experiment. Two points have to be addressed: i) the contributions of El retardation 
and higher multipoles, which are not contained in S^^^, but which will contribute to 
the experimental cross section and ii) the contribution of the high energy tail. 
Regarding point i), additional effects of the El retardation and higher multipoles have 
been calculated in Ref. (Golak et al., 2002) for the AV18+UIX potentials in a Faddeev 
calculation. Using those results one finds that El retardation and higher multipoles 
increase ?7i_i(135) by about 1% only. In fact according to Gerasimov(Gerasimov, 
1964) there is a large cancelation of El retardation and other multipoles. There is no 
reason for a larger effect in ^He. On the contrary, since the leading isovector magnetic 
dipole (Ml) transition is suppressed in this nucleus (it is zero for an S-wave) one can 
expect an even smaller contribution. As to point ii), considering the fall-off of the 
triton ElUR cross section around pion threshold from Ref. (Golak et al., 2002) one 
obtains a result very similar to our ^He case, namely p ^ 1.5. However, including the 
other multipole contributions one gets a considerably smaller value, namely p ~ 1.1. 
This is no contradiction with the small increase of 1% for S^^^, since the inverse 
energy weighted cross section integrated from 100 to 135 MeV gives only a rather 
small contribution to the sum rule. On the other hand one would overestimate E^^^ 
taking the full cross section. Thus it is suggested that the tail contribution to an 
experimental BSR should be estimated using the theoretically established fall-off co'"^ 
with p ^ 1.5. 



n 



Neutrino Scattering on 
Light Nuclei in Supernovae 

The inelastic scattering of neutrino off A = 3 nuclei and ^He is calculated microscop- 
ically at energies typical for core collapse supernova environment. The calculation 
is carried out with the Argonne fig nucleon-nucleon potential and the Urbana IX 
three nucleon force. Full final state interaction is included via the Lorentz integral 
transform (LIT) method. The contribution of axial meson exchange currents to the 
cross sections is taken into account from effective field theory of nucleons and pions 
to order 0{Cf'\ The main results of the chapter are based on a series of papers 
pubhshed in Gazit and Barnea (2004, 2007b,a,c); O'Connor et al. (2007). 

The chapter starts with an overview on core collapse supernovae, emphasizing the 
abundances and role of ^He and A=3 nuclei in the processes taking place in the SN 
environment. The cross-sections are then presented and discussed. Finally, neutrino 
mean-free-paths in the neutrinosphere region are evaluated. 

7.1 Core-collapse Supernovae 

A massive star, above about 10 solar masses, evolves over millions of years, to create 
an onion-like structure of layers, with lighter elements in the atmosphere and heavier 
in the center, where the gravitational pressure is highest. The cycle of fusion stops 
when the core burns Si nuclei to iron. Further burning can not be a source of energy, 
as the binding energy per nucleon is peaked in iron. Consequently, the evolution of 
the star comes to an end. Observations have shown that this death is not peaceful, 
it is followed by an enormous explosion, shining in a magnitude of a galaxy, believed 
to occur in the following chain of events. 

When the core reaches the Chandrasekhar mass it cannot support its own gravity 
and collapses. The collapse halts only when the core reaches nuclear density, in which 
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short-range nuclear forces dominate the dynamics. This initiates an outgoing shock 
wave. As the shock travels through the outer layers of the core it losses energy due 
to thermal neutrino radiation, and dissociation of iron nuclei (bound by about 8 
MeV/nucleon) to a particles and free nucleons. 

The fate of this shock is a matter of debate. In the "prompt shock" mechanism, 
the shock has enough initial energy to burst through the core and cause an explosion. 
However, in modern calculations, about 100 miliseconds after bounce the shock stalls 
at about 200 km from the center and becomes an accretion shock, with matter from 
the outer layers of the star falling through it. During this period the newly born 
proto-neutron star in the center of the star cools by neutrino emission. This huge 
burst of neutrinos, in total energy of almost 10^^ ergs, is believed to deposit enough 
energy in the matter below the shock to revive the shock. This so called "delayed 
shock" mechanism has not been quantitatively proved in numerical calculations, but 
received a qualitative proof in the observation of SN 1987A. Comparing the energy 
of the neutrino-burst to the gravitational well of the star, it is easy to realize that 
one needs only 1% of neutrino energy to be transferred to matter. This turns to be 
a highly delicate problem, in which the accuracy of the microscopic input is a crucial 
factor. The heating is achieved through elastic scattering on nucleons and electrons. 
Haxton (1988) has suggested that inelastic scattering on nuclei can increase neutrino- 
matter coupling. The higher energy of the heavy flavored neutrinos^ favors neutral 
scattering on nuclei as an energy transfer method. Thus, neutrino interaction with 
nuclei, which exist below the shock, is believed to have an important part in the 
explosion process. This will be discussed thoroughly in Subsec. 7.1.1. 

Another outcome of the huge neutrino flux is changing the chemical evolution of 
the star, mainly by break-up of nuclei into fragments of nuclei which act as seeds 
for different nucleosynthesis processes. ^He break-up by neutrinos has a role in this 
process, which will be reviewed in Subsec. 7.1.2. 



^The characteristic temperatures of the emitted neutrinos are about 6 — 10 MeV for i/^ ,- (t'p,T), 
5 - 8 MeV for i>p, and 3 - 5 McV for i^p 
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7.1.1 Shock revival and Light Nuclei 

In order to estimate how significant is neutrino interaction with hght nuclei in the 
shocked region, one needs to know not only the neutrino luminosity but also the 
mean-free-paths of neutrinos in the region. Thus one needs the composition of the 
shocked regions and the cross-sections for neutrino scattering on light nuclei. In 
Fig. 7.1 one can view the density profile of an IIMq progenitor, 120 miliseconds after 
bounce (Livne, 2007). The calculation is ID simulation, with full neutrino-transport, 
and using the equation of state developed by Shen et al. (1998). The red line in 
the Figure is the a mass-fraction, which is the most abundant nuclei in the shocked 
region by this calculation. The rest of the composition are nuclei and free electrons. 
One can see that in this simulation the amount of ^He is about 1% up to 100 km off 
center and rises to around 10 — 20% just below the shock. Thus neutrino scattering off 
'^He could play a role in the revival of the shock, as first indicated by Haxton (1988). 
Moreover, the energy transfer due to elastic scattering is low, lo ~ T"^ /m, therefore 
inelastic scattering could be as important as elastic scattering on nucleons. One has 
to keep in mind that with every inelastic scattering there is a substantial release of 
energy - larger than the break-up energy of 20 MeV. A complete study on the role 
of ^He excitation on the shock revival was not done, and is called for. However, a 
first study was carried out by (Ohnishi et al., 2006), and showed a small effect on 
accretion shock instabilities. 

Nuclear statistical equilibrium (NSE) models, like those used by Shen et al. 
(1998), predict abundances based on binding energies and the quantum numbers 
of nuclei. However, NSE models only treat approximately (or neglect) strong in- 
teractions between nuclei, and consequently break down as the density increases. 
Recently a description of low-density nuclear matter (composed of neutrons, pro- 
tons and alpha particles) in thermal equilibrium based on the virial expansion was 
developed(Horowitz and Schwenk, 2006a,b). The virial equation of state (EoS) sys- 
tematically takes into account contributions from bound nuclei and the scattering 
continuum, and thus provides a framework to include strong-interaction corrections 
to NSE models. This formalism makes model-independent predictions for the con- 
ditions near the neutrinosphere, i.e. for densities of p ~ 10^^~^^g/cm^ and high 
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distance from center [km] 



Figure 7.1: Density profile of IIMq star, 120 miliseconds after bounce. The black 
line is the total density, and the red line is the a mass fraction. One can view the 
accretion shock at about 200 km off center. 



temperatures T ~ 4MeV (Costantini et al., 2004; Lunardini and Smirnov, 2004). 
In particular, the resulting alpha particle concentration differs from all equations of 
state currently used in SN simulations, and the existence of trinuclei in this area is 
expected. For example, in Fig. 7.1 the region of about 10 — 30 km off center is char- 
acterized by this density range. Thus, the calculation substantially underestimates 
the abundances of A = 3, 4 nuclei in this region. 

To determine the abundance of A = 3 nuclei near the neutrinosphere in super- 
novae, we presented in O'Connor et al. (2007) an extension of the virial approach 
to the EoS to explicitly include neutrons, protons, a particles, ^H and '^He nuclei; 
deuterons are included as a bound state contribution to the proton-neutron virial 
coefficient. 

The resulting mass fractions Xi = Airti/ni^ {jii is the density of the i element of mass 
Ai, and Ub is the baryon density), are shown in Fig. 7.2 for neutrinosphere densities 



and temperatures, and various proton fractions Yp 
the ratio of proton number to baryon number. 



[nr 



2na + 2n3He + n3ii)/nb, i.e. 
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Figure 7.2: Mass fractions of nucleons and A = 3, A nuclei in chemical equilibrium 
as a function of temperature T. The top and bottom rows correspond to a density 
of 10^^ g/cni'^ and 10^^ g/cm^ respectively, and from left to right the proton fractions 



are Yp = 0.1, 0.3 and 0.5. 



The main result of this EoS is the model-independent prediction that the mass- 
three fraction can be significant (up to 10%) near the neutrinosphere. Thus, we 
present in Sec. 7.3.4 predictions for the neutral-current inclusive inelastic cross-sections 
on mass-three nuclei, based on microscopic two- and three-nucleon interactions and 
meson-exchange currents, including full final-state interactions via the Lorentz inte- 
gral transform (LIT) method. 



7.1.2 Nucleosynthesis and u — a Interaction 

Due to the high temperature of /i and r neutrinos (and anti- neutrinos), a substan- 
tial amount of them carry more than 20 MeV, and may dissociate the ^He nucleus 
through inelastic neutral current reactions. This creates the seed to light element 
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nucleosynthesis in the surrounding stellar envelope (Woosley et al., 1990). As a part 
of this explosive nucleosynthesis, or "z/ -nucleosynthesis" , a knock out of a nucleon 
from a ^He nucleus in the helium rich layer, followed by a fusion of the remaining 
trinucleus with another a particle, will result in a 7-body nucleus. This process is 
an important source of ^Li, and of ^^B and ^^F through additional a capture reac- 
tions. In fact the production of ^^B in supernovae (also due to neutrino interaction on 
^^C) can be calibrated by measuring ratio of ^^B to ^°B in meteorites. Woosley and 
Weaver (1995) have suggested that the calculated abundance is too large. Yoshida 
et al. (2005) have shown that the Li and B production depends highly on the neutrino 
temperature, and mainly on the spectrum character of the high energy tail (due to 
the high dissociation energy of a). They used this observation to fit the neutrino 
temperature to reproduce the measured abundance. However, a correct description 
of the process must contain an exact, energy dependent cross-section for the neutral 
inelastic a — v reaction, which initiates the process. 

The relatively low temperature of v^^ and v^ emitted from the core suppress the 
probability for inelastic reactions of these neutrinos with ^He in the supernova sce- 
nario. However, the mixing of the first and the third neutrino mass eigen-states is 
enhanced, due to matter effects, in the Oxygen/Carbon layers, just below the He 
layer. This significantly changes the spectrum of electron (anti-) neutrinos reaching 
the He layer, and yields a secondary source of energetic electron neutrinos. As in- 
dicated by Yoshida et al. (2006), the resulting charge current inelastic reactions on 
^He could double the aforementioned yields. In this reference, they speculate that 
a test for neutrino oscillation properties could be observing the ratio of ^Li to ^^B 
in stars. Yet again, this could lead to quantitative conclusions only if neutrino-^He 
cross-sections are available and accurate. 

Elements heavier than iron are produced mainly by a different nucleosynthesis 
process - neutron capture. The so called r-process, in which this capture is rapid, 
is still missing a neutron-source in nature. Epstein et al. (1988) have suggested that 
neutrons produced by neutrino spallation of ^He in the helium rich layer could be 
such a source. 

However, the most probable source of r-process neutrons is the neutrino-driven 
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wind blowing from the surface of a newly born proto-neutron star (PNS). This is a 
high entropy region, thus all protons combine with neutrons to form ^He, with excess 
neutrons left over. This makes the r-process picture puzzling, as the free neutrons 
that initially accompany the ^He, are then converted to protons by charge current 
reactions, which then form more ^He (Woosley and Hoffman, 1992). Thus it is very 
hard to maintain a neutron excess (necessary for the r-process) in the wind. Neutrino 
cross sections for ^He breakup are part of the detailed balance that determines what 
residual neutron abundance will remain in a neutrino-driven wind. 

7.2 Weak Interaction Observable — Beta Decay Rate 

In order to give parameter free predictions of weak processes, one needs the LECs of 
the EFT. There is only one unknown LEG, which is calibrated to fit a weak interaction 
observable. 

A most researched weak interaction observable is the (3 decay. In this decay, a 
neutron (which can be free or bound in a nucleus) spontaneously transforms into a 
proton, accompanied by electron and anti-neutrino emission. This instability could be 
found in many elements, transforming one element to another. The lightest elements 
undergoing this decay are: n — > p+e~ + z/e, ^H ^ B.e+e~ + Ue and ^He — ^^ Li+e~ + Ue. 

The /5 decay rate provides a direct test for the weak interaction model and the 
axial currents in the nucleus. This aspect is even more appealing as the /5 transitions 
are between bound-states of nuclei, and usually include small momentum transfer. 
Thus, the decay is governed by the Gamow-Teller and Fermi operators, which are 
proportional to the low momentum transfer approximation of E^ and Cq operators. 
In fact, to leading order the f3 decay rate is simply given by (see for example (Schiavilla 
and Wiringa, 2002)) 



G'(+) 



2^5 



r?. = /.^^^(|Fr + ^l|GTn (7.1) 
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where the Gainow Teller and Fermi reduced matrix elements: 

A 

i=l 
A 

GT = {f\\Y.^^,+^\^. (7.3) 

ft includes the entirety of the kinematics: 

j-WH-rnJ 

ft= deeV^^^{m-mJ-efF^+\Zf,e) (7.4) 



here e is electron energy in units of electron mass, m are masses in the same units, 
with subscript i{f) indicating initial (final) state, respectively. Since isospin mixing is 
usually small for light nuclei, the Fermi operator can be evaluated accurately. Thus, 
by using this formula, one can deduce the experimental GT to very good accuracy 
(Chou et al., 1993): 

GT{n-^p + e- + Ue) = V3-(l± 0.003) (7.5) 

GT(3H ^^ He + e" + z/e) = Vs ■ (0.961 ± 0.005) (7.6) 

The axial coupling constant, qa = 1.2670±0.0004, is calibrated using the neutron life- 
time, where the matrix element is analytic. For the triton, axial MEG in the nuclei 
influence the half-life. Thus, the triton GT is used to calibrate unknown constants 
in the MEG model. In the case of MEG based on the EFT introduced in Ghap. 3, 
the only unknown is the renormalization LEG dj.. 

Keeping the philosophy of the EFT* approach, the AV18/UIX Hamiltonian is 
used as the nuclear interaction for the calculation of the trinuclei wave functions. 
One then calculates the GT {E\) using the EFT based MEG, and calibrates dr to 
reproduce the half life. This procedure follows (Park et al., 2003). One could use 
their result for dr{A). However, as a test for both the nuclear wave functions and the 
numerics. 

When reproducing, there are additional benchmarks one should note. First, the 
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Method 


Binding Energy [MeV] 


^H 


^He 


EIHH 

CHH 

FY 


8.471(2) 
8.474 
8.470 


7.738(2) 
7.742 
7.738 


Experimental 


8.482 


7.718 



Table 7.1: Binding energies of ^H and ^He calculated using AV18/UIX Hamiltonian 
model compared to the same calculation done by using FY equations and expansion 
on correlated hyperspherical harmonics (CHH) basis (Nogga et al., 2003). For the 
EIHH calculation, the number in parenthesis indicates the numerical error. Also 
shown are the experimental values. 

calculation can be checked by comparing binding energies of the trinuclei to those 
achieved by the FY and the correlated hyperspherical harmonics methods (Nogga 
et al., 2003), with the same potential model. This appears in Table 7.1. One can see 
a very good comparison between the calculations. 

An additional benchmark is a calculation of the single operator contribution to 
the GT amplitude of ^H half life (which will be denoted by GTi^) (Schiavilla et al., 
1998) with the same potential model, using the CHH method. This comparison reads: 
GTib(CHH) = 1.596 ^ GTib(EIHH) = 1.598(1). 

Finally, I reproduce the dr{A) of (Park et al., 2003). My results are: 



dr = 1.02 ± 0.02 ± 0.08 for A = 500 MeV 

dr = 1.78 ± 0.03 ± 0.08 for A = 600 MeV 

dr = 3.8 ± 0.06 ± 0.12 for A = 800 MeV 



(7.7) 



Where the first error is numerical, i.e. convergence of the calculation, and the second 
correspond to the experimental uncertainty in the triton GT strength. This is to be 
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compared with (Park et al., 2003), 

dr = 1.00 ± 0.07 for A = 500 MeV 
dr = 1.78 ± 0.08 for A = 600 MeV 
dr = 3.9 ± 0.10 for A = 800 MeV. 

The minor differences are of the order of those encountered in other observables. 
The differences in the error-bars are probably due to differences in the experimental 
source. 

Summing up, the results of this section give us the needed information to carry 
on parameter free predictions of weak processes. As a side benefit, I get a verification 
of the numerical methods used, and the derived weak currents. 

7.3 Neutrino Scattering Cross-sections Calculation 

The calculation of the neutral-current inclusive inelastic cross-sections follows the 
procedure used to calculate the photoabsorption process. We solve the bound state 
problem based on the Argonne fig nucleon-nucleon potential (Wiringa et al., 1995) 
and the Urbana IX three- nucleon interaction (Pudliner et al., 1997). The bound state 
binding energies of the trinuclei are given in the previous section, and those of ^He 
in Table 6.1. 

Neutrino scattering on A = 3 and ^He nuclei only induces transitions to continuum 
states, since the nuclei have no bound excited states. The correct description of the 
process is achieved via the LIT method. The resulting Schrodinger-like equations are 
solved using the effective interaction hyperspherical harmonics (EIHH) approach of 
Chap. 5. 

We use chiral EFT meson-exchange currents (MEC) as presented in Chap. 3. The 
EFT approach is an appropriate approximation, as Q ^ 60 MeV is the typical energy 
in our processes of interest, and the cutoff A is of the order of the EFT breakdown 
scale A = 400 — 800 MeV. In configuration space, the MEC are obtained from a 
Fourier transform of propagators with a cutoff A (cf. Eq. (3.36)). This leads to a 
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cutoff dependence, which is renormahzed by the cutoff-dependent counterterm dr{A), 
given in Eq. (7.7). 

7.3.1 Leading Multipole Contributions 

In the supernova scenario one has to consider neutrinos with up to about 60 MeV. The 
main resuh of this is that it is sufficient to retain contributions up to 0{q^) in the mul- 
tipole expansion. This conclusion is checked on ^He, and is accurate to the percentage 
level. This is consistent with the discussion of Sec. 2.3. The leading contributions to 
the inelastic cross-section are due to the axial vector operators E^, E^, M^, L^, Lq 
and the vector operators C^ , EY , L\ . As explained in the previous chapter, the lead- 
ing vector operators are all proportional to each other due to the Siegert theorem, 
valid in this energy range (see discussion in Sec. 2.3). The axial charge operator Cq 
and the magnetic M^ operator are both proportional to the inverse of the nucleon 
mass - thus have no substantial contribution, exactly as relativistic corrections. In 
fact, for a percentage level accuracy, the following long wavelength approximation can 
be used (Donnelly and Walecka, 1976), for the impulse approximation the spin-orbital 
parts of the operators are 



C^id) 



A-K 



.Ar ^ ^ 



v47r 



CUq) = y>^iM(r) 
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2 



E2KM) = M/ 5^a|[c? ® Y,if)]^, 



L2Mi<l) = y ^KAq) 

The axial MEC operator, given in Eq. (3.48), is invariant to an exchange between 
the particles. Thus, it does not influence multipoles of opposite symmetry. The MEC 
will have a substantial influence only on the GT {E\) operator, as it will be discussed 
in the following subsections. 

7.3.2 Inelastic Cross-sections and Energy-transfer 

The result of the current calculation gives an energy and angle dependent inclusive 
cross-section. In order to present the cross-sections, we use the fact that the neutrino 
spectra are approximately thermal. Thus, the calculated cross sections are averaged 
over energy and angle, assuming a Fermi-Dirac distribution for the neutrinos with 
zero chemical potential, temperature T^, and neutrino momentum A;, 

where A^~^ = 2 ^^j^(— l)'""*"^/n^ is a normalization factor. The quantities of inter- 
est are the cross-sections and energy transfer cross-sections averaged over neutrino 
spectra of temperature T^, 

poo /• 1 

{a)n= / dio dhf{n,h) — , (7.9) 

/•CxD p 1 

{uja)T, = / duj dkif{T^,ki)uj——, (7.10) 

where kij are the initial and final neutrino energy, u = ki — kf is the energy transfer, 
and tUth denotes the threshold energy of the breakup reaction. For ^He, ^He and ^H 
the threshold energies are Utt = 19.95, 6.23, 5.5 MeV respectively. 
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It is well known that realistic 2-body NN potentials lead to an under-binding of 
about 0.5—1 MeV for the ^He and the triton nuclei and an under-binding of about 3—4 
MeV for ^He. The AV18 threshold energies are 16.66 MeV, 5.37 MeV and 4.67 MeV 
for ^He, ^He and triton breakup reaction, respectively. Thus the AV18 model has a 
discrepancy, A ^ 3.3 MeV for the threshold energy of the ^He, A ^ 0.83 MeV for 
^He break-up, and A ^ 0.86 MeV for '^H, with respect to the experimental inelastic 
reaction threshold. 

Thermal averaging is highly influenced by the threshold energy, thus very sensitive 
to this discrepancy. In order to correct for this difference the response function is 
shifted to the true threshold, i.e. R{uj) — > Riu — A), when only NN potentials are 
used. 

7.3.3 Neutrino scattering on ^He nucleus 

In terms of "many-body" nuclear physics, '^He is a closed shell nucleus, i.e. its total 
angular momentum is zero and its spin structure is governed by an S-wave. Some 
of the leading operators are strongly suppressed due to the special structure of the 
ground state. 

In fact, the Gamow- Teller operator contributes only due to the small P and D 
wave components of the ground state wave function. This reflects on the axial MEG 
contribution to the cross-section. In magnitude, the Gamow-Teller strength is less 
than 0.25% of the total cross-section, for all considered nuclear potentials: AV18, 
AV8', and AV18/UIX. Due to MEG, the strength grows to 1% ± 0.5%. The error is 
due to the A dependence. In this case, the part of the error which correspond to the 
triton half-life is of no importance. The MEG contribution is calculated only for the 
AV18/UIX potential model, and dr is taken from Eq. (7.7). 

The special internal structure of the nucleus affects also the MY operator, which 
vanishes for ^He. In addition, ^He is an almost pure zero-isospin state (Nogga et al., 
2002), hence the Fermi operator vanishes. 

As a result, the leading contributions to the inelastic cross-section are due to the 
axial vector operators i?^, M^, L^, Lq and the vector operators C^ ^ EY , L^ , which 
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Figure 7.3: Relative error in the sum-rule of the leading response functions with 
respect to the hyper-angular momentum quantum number K (for ^He). Calculated 
using the AV8' NN potential model. The error bars reflect the uncertainty in inverting 
the LIT. 



are all proportional to the momentum transfer q. 

The combination of the EIHH and LIT methods brings to a rapid convergence in 
the Response functions. In Fig. 7.3, one can see the relative error in the sum-rule of 
the main response functions when calculated with NN potential (in this case Argonne 
Wg, but the conclusions hold for AV18 as well), with respect to the hyper-angular 
momentum quantum number K. It can be seen that upon convergence the relative 
error is well below 1%. The error bars presented reflect the error in inverting the LIT. 
This convergence pattern holds also for the AV18 potential model. Bearing in mind 
that the cross-section, up to kinematical factors, is the sum of the response functions, 
this is a measure of the accuracy in the calculated cross- sect ion. 

The convergence rate is outstanding even when including the 3NF. In Fig. 7.4 
we present for the leading multipoles the convergence of the LIT as a function of 
K. It can be seen that the EIHH method results in a rapid convergence of the LIT 
calculation to a sub-percentage numerical accuracy. We conclude that the 3NF does 
not affect the convergence rate of these operators. 
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Figure 7.4: Convergence of ^o^o^/l for the leading operators, as a function of the 
HH grand angular momenta K (for ^He). 



The general form of the cross-section is shown in Fig. 7.5. 

In Table 7.2 we present the temperature averaged total neutral current inelastic 
cross-section as a function of the neutrino temperature for the AV8', AV18, and the 
AV18+UIX nuclear Hamiltonians and for the AV18+UIX Hamiltonian adding the 
axial MEC. From the table it can be seen that the low-energy cross-section is rather 
sensitive to details of the nuclear force model. This sensitivity is visible already when 
comparing the AV8' and AV18 results, where an effect of about 10% is discovered. 

The 3NF have a bigger influence of about 25%. Similar tendency was also ob- 
served for the hep process (Marcucci et al., 2000). The effect reduces with neutrino 
temperature, from 35% for % = 2 MeV, to 23% for % = 12 MeV. 

The results are of the same order of magnitude as previous estimates (Woosley 
et al., 1990), also appearing in the Table, though the differences can reach 25%. The 
current work predicts a stronger temperature dependence, with substantial increment 
at high temperatures. This indicates a different structure of the predicted resonances. 

The cross-section is dominated by the axial E2 and Mi multipoles, which for 
example exhaust about 90% of the cross-section at Tj^ = 10 MeV. 
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T 




(-^)r = ii« + 


O^ 10-^2^^. 




MeV] 


AV8' 


AVI 8 


AV18+UIX 


AV18+UIX+MEC 


Woosley (1990) 


4 


2.09x10-3 


2.31x10-3 


1.63x10-3 


1.66x10-3 


— 


6 


3.84x10-2 


4.30x10-2 


3.17x10-2 


3.20x10-2 


3.87x10-2 


8 


2.25x10-1 


2.52x10-1 


1.91x10-1 


1.92x10-1 


2.14x10-1 


10 


7.85x10-1 


8.81x10-1 


6.77x10-1 


6.82x10-1 


6.87x10-1 


12 


2.05 


2.29 


1.79 


1.80 


1.63 


14 


4.45 


4.53 


3.91 


3.93 


— 



Table 7.2: Temperature averaged neutral current inclusive inelastic cross-section per 
nucleon as a function of neutrino temperature. The last column is the calculation in 
Woosley et al. (1990). 



T 


{a)T [10-^2^^2] 


[MeV] 


i^x^V^) 


{^.,K) 


(^^6,6 ) 


(^^6,6+) 


2 
4 
6 
8 
10 


1.47x10-^ 
1.73x10-3 
3.34x10-2 
2.00x10-1 
7.09x10-1 


1.36x10-^ 
1.59x10-3 
3.07x10-2 
1.84x10-1 
6.54x10-1 


7.40x10-^ 
8.60x10-3 
1.63x10-1 
9.61x10-1 
3.36 


5.98xl0-« 
6.84x10-3 
1.30x10-1 
7.68x10-1 
2.71 


T 


{auj)T [10-42 MeVcm2] 


[MeV] 


(^.,^i) 


{^xK) 


{j^e,e ) 


(^e,e+) 


2 
4 
6 
8 
10 


3.49x10-5 
4.50x10-2 
9.26x10-1 

5.85 
21.7 


3.23x10-5 
4.15x10-2 
8.56x10-1 

5.43 

20.2 


1.76x10-4 
2.27x10-1 

4.56 

28.4 

103.8 


1.42x10-4 
1.80x10-1 

3.70 

22.9 

84.4 



Table 7.3: Temperature averaged inclusive inelastic cross-section (upper part) and 
energy transfer cross-section (lower part) per nucleon as a function of temperature. 
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Figure 7.5: Temperature averaged inelastic cross-sections at temperature T = 10 
MeV for ^He. The solid line is the differential cross-section, (^)t = \\{^ + S)t, 
(left scale). The dashed line is the differential energy transfer cross-section, (i^^)t = 



1 j_ 

2 A 



{lo 



t7 + ^S)r'("ght scale). 



In Table 7.3 we present (for AV18+UIX+MEC) the temperature averaged cross- 
section and energy transfer as a function of the neutrino temperature for the various 
processes. In both tables it can be seen that the charged current process is roughly a 
factor of five more efficient than the neutral current process. 



7.3.4 Neutrino scattering on A = ?> nuclei 

Contrary to ^He, the GT is not suppressed for the trinuclei, which have J = \ total an- 
gular momentum. At low-momentum transfer, the Gamow- Teller operator dominates 
for the cross section. It contributes about 60% of the cross-section at a temperature 
of 1 MeV. At higher- momentum transfer, higher-order multipoles (mainly the axial 
E2 and Ml) start to play an important role, thus the GT part of the cross-section 
decreases to about 20% at Tj, = 5 MeV, and to about 10% at neutrino temperature 
of 10 MeV. 

As a consequence of the descending importance of GT, the MEC contribution 
also decreases with temperature. In Fig. 7.6, one can view the relative contribution 
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n [MeV] 


^H 


3He 


1 


1.97x10-6 1.68x10-5 


3.49x10-*^ 


2.76x10-5 


2 


4.62x10-^ 4.73x10-3 


6.15x10-^ 


5.94x10-3 


3 


5.53x10-3 6.38x10-2 


6.77x10-3 


7.41x10-2 


4 


2.68x10-2 3.37x10-1 


3.14x10-2 


3.77x10-1 


5 


8.48x10-2 1.14 


9.70x10-2 


1.25 


6 


2.09x10-1 2.99 


2.35x10-1 


3.21 


7 


4.38x10-1 6.61 


4.87x10-1 


7.03 


8 


8.20x10-1 13.0 


9.03x10-1 


13.7 


9 


1.41 23.4 


1.54 


24.6 


10 


2.27 39.3 


2.47 


41.2 



Table 7.4: Averaged neutrino- and anti-neutrino-3H and -3He neutral-current inclu- 
sive inelastic cross-sections per nucleon (A=3), {cr)T^ = jj ( o-^ + a-ij)^^ (left columns), 
and energy transfer cross-sections, {uja)T^ = j^ {uja^ + uja-ij)^^ (right columns), as a 
function of neutrino temperature Ty, in units of 10-^2 cm^ and 10-^2 MeVcm2 respec- 
tively. 



of MEC to the energy transfer, 






. The width of the line indicates the 



error in the theoretical estimate, due to the cutoff dependance. The cutoff dependence 
of both cross-sections is ~ 2% for 1 MeV and < 1% for higher temperatures. This is 
a rather strict validation of the calculation. The numerical convergence is very good 
for the trinuclei ~ 1%. 

The 3NF affect the cross-sections by about 15%, which is less than the their effect 
in ^He, however still substantial. 

While not directly important for the shock revival, the asymmetry between the 
scattering of neutrinos and anti-neutrinos increases with temperature: the difference 
in the energy transfer grows gradually from 3% for a neutrino temperature of 3 MeV 
to > 50% for 10 MeV temperatures. 



7.3.5 Neutrino energy loss due to inelastic scattering 

We can combine the energy transfer cross-sections with the A = 3, 4 mass fractions of 
Fig. 7.2 to calculate the neutrino energy loss due to inelastic excitations of A = 3, 4 
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Figure 7.6: 

{ua)^° MEC 



Relative contribution of MEC to the energy transfer, i.e. 
The width of the hne indicates the error in the theoretical es- 



timate, due to the cutoff dependance. 



nuclei. In Fig. 7.7 we compare the total energy transfer cross-sections for ^H, ^He and 
^He. The comparison clearly shows the huge difference between the cross-sections, 
originating in the "open-shell" character of the trinuclei. 

The neutrino of energy E^ will decrease due to inelastic excitations, and heat the 
matter, at a rate dE^/dx given by 



dE, 
dx 



i=3H, 3He, *Hc 



(7.11^ 



For simplicity, we neglect the energy transfer from nuclei to neutrinos required by 
detailed balance. This is strictly correct only in the limit T ^ Tj,. 

In Fig. 7.8, the neutrino energy loss due to inelastic scattering is shown for a 
density of 10^^ g/cm^ and neutrino temperature T,^ = 6 MeV, as a function of the 
matter temperature for various proton fractions. For T > 4 MeV, the energy loss is 
dominated by the contributions from ^H nuclei. The total abundance of A = 3 nuclei 
depends only weakly on the proton fraction (see Fig. 7.2), which is reflected in the 
weak dependence of the neutrino energy loss as a function of proton fraction. Finally, 
for lower densities, mass-three nuclei are less abundant (see Fig. 7.2), and therefore 
also their contributions to the neutrino energy loss. 
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Figure 7.7: Neutrino total energy transfer cross-section a function of the neutrino 
temperature T^, for A = 3, 4 nuclei. 
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Figure 7.8: Neutrino energy loss dE^jdx for inelastic excitations of A = 3,4 nuclei 
as a function of the matter temperature T at a density of lO^^g/cm^. We assume 
that the neutrino energies are characterized by a Fermi-Dirac distribution with a 
temperature T,^ = 6 MeV. The contributions from ^H, ^He and ^He nuclei, and the 
total neutrino energy loss are shown for proton fractions 1^ = 0.1, 0.3 and 0.5. 




Conclusions 



In this thesis I have presented the methods and the resuhs of ab-initio calculations 
of electro-weak reactions on light nuclei. The reactions are calculated using realistic 
NN potential, Argonne wig, and includes the Urbana IX 3NF. The full interaction is 
taken into account not only for the ground state, but also for the continuum states 
via the LIT method. For the solutions of the differential equations I use expansions 
in hyperspherical harmonics via the EIHH approach. 

In the next two sections I will briefly summarize the main results and implications 
of the work. 

8.1 Photoabsorption on ^He 

I have presented the first complete calculation of the ^He total photoabsorption cross 
section using a realistic nuclear force (AV18 NN potential and the UIX-3NF). The 
results show a rather pronounced giant dipole peak typical to many body nuclei. This 
microscopic calculation of a complex nucleus settles the long living debate regarding 
the character of the peak. 

The effect of the 3NF reduces the peak height by only 6%, less than expected 
considering its large effect of almost 20% on the ^He binding energy and its different 
role in the three-nucleon system. Beyond the giant dipole resonance 3NF effects 
become much larger. With growing energy transfer their importance increases and at 
pion threshold one finds a cross-section enhancement of 35%, about twice the effect 
discovered in ^H/'^He photoabsorption. 
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The calculated cross-section was integrated with energy-transfer weights to inves- 
tigate sum-rules for ^He, which can be connected to important ground state proper- 
ties. Three well known photonuclear sum rules for ^He are studied, viz the Thomas- 
Reiche-Kuhn, the bremsstrahlung and the polarizability SR. Two new equivalences 
for the BSR have allowed deduction of information about two-body properties of the 
nuclear ground state, like the proton-proton, neutron-neutron and proton-neutron 
distances. 

A result with some beauty is tying the SR to the configuration symmetry of the 
nucleus. In particular, the configuration tetrahedral symmetry of ^He was tested, 
and found this symmetry to be slightly broken. An experimental way to access 
this symmetry breaking is proposed, via a measurement of the BSR which could be 
performed in one of the existing or planned low-intermediate energy photonuclear 
facility. 

The calculated cross-section can be compared to numerous experiments estimat- 
ing the reaction. Close to threshold the theoretical cross section agrees quite well 
with experimental data. In the giant resonance region, where there is no established 
experimental cross section, the results are in good agreement with the data of Nilsson 
et al. (2005) and Wells et al. (1992), while a strong disagreement is found when com- 
pared with the data of Shima et al. (2005). In order to understand whether a nuclear 
force model, which is constructed in the two- and three-nucleon systems, is sufficient 
to explain the four-nucleon photodisintegration, further experimental investigations 
are mandatory. 

8.2 Neutrino Interaction with ^He and A=3 nuclei 
in Supernova 

The inelastic scatterings of neutrino off ^He, ^H and ^He are calculated microscopi- 
cally at energies typical for core collapse supernova environment. These are the first 
microscopic estimates of the processes, carried out with realistic nucleon forces. All 
breakup channels and full final-state interactions were included via the LIT method. 
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The contribution of axial meson exchange currents to the cross sections is taken into 
account from effective field theory of nucleons and pions to order 0{Q^), which re- 
produce the triton half-life. 

I estimate the overall accuracy of the calculation to be of the order of 5%. This 
error is mainly due to the strong sensitivity of the cross-section to the nuclear model, 
in particular to the 3NF. The numerical accuracy of our calculations is of the order 
of 1%. The contribution of the axial MEC for ^He is small, in the percentage level. 
However, for the trinuclei the MEC contribution could be up to 30% for low tem- 
perature. The cutoff dependence for both A = 3, 4 nuclei, which is indicative to the 
validity of the hybrid EFT approach, is negligible, less than a percent for neutrino 
temperature above 2 MeV. 

Using the virial abundances and the microscopic energy transfer cross-sections, 
it is found that mass-three nuclei contribute significantly to the neutrino energy loss 
due to inelastic excitations for T > 4 MeV, conditions typical to the neutrinosphere. 
The predicted energy transfer cross-sections on mass-three nuclei in this area are 
approximately one order of magnitude larger compared to inelastic excitations of ^He 
nuclei. 

To fully assess the role of neutrino breakup of ^He in core-collapse, the cross- 
sections should be included in SN simulations. The u — a interaction is of central 
importance also for correct evaluation of the z/— and r— nucleosynthesis processes. 

With the present work, I have made an important step in the path towards a more 
robust and reliable description of neutrino-nuclei interaction role in core-collapse 
supernovae. 




Weak Currents In the 



Standard Model 



In this appendix, the general structure of the weak currents in nuclei are deduced 
from the standard model. 

In the case of neutral neutrino scattering off a point proton, making use of the 
Feynman diagrams, with the Z^ boson propagator approximated as g^j^/M"^ one gets: 



Tf, = -M^'-u{kj) 



-w 



AcosO 



w 



-r{l-l5) 



u{ki) 



(A.l) 



-w 



Acosdw 



7^((l-4sin2^H^)-75) 



■u{Pf) 



u{pi 



u{p2)l^i ( -(1 - 75) -2siv?ew ] u{pi] 



Where G = „. ,a ^^„ is the Fermi constant. A similar calculation for neutrino neutral 
scattering on neutron yields a different result: 



Tf, = ^[u{h)r{l - l,Hh)] 



'"(^2)7/. ( -2(1-75) ) u{pi) 



(A.2) 



One can combine the results of equations A.2 and A.l, using the isospin symmetry 
to get for the nucleon current: 



J'' = U{P2 



2 ' ^ 



75)-2sin2^H^7'^-(l + ro) 



u{pi 



(A.3) 



For a complex nucleon, one keeps the symmetries of this current to write 



7(0) 



(1 - 2 . sin^ 9w)^Ji: + ^J^ - 2 . sin^ ^^1 jj 



(A.4) 
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Where J^ and Jrf have vector and axial symmetry, respectively. The second equality 
takes advantage of the conserved vector current hypothesis, by which the weak-vector 
current is an isospin rotation of the electro-magnetic current. 
A similar calculation for the charged current yields: 

J^^' = ^4iJ: + J^) (A.5) 




Lepton Current 



In this appendix the calculation of the lepton current matrix element is demonstrated. 
I will show this for the case of a neutrino scattering, as calculation of this matrix 
element for /5 decay or electron capture follow the same line. 
The leptonic current is, 

j,{x) = Mx)j,{l + 75)V^.'(x) = l^e-'"^-'' (B.l) 

Where, for lepton reaction 

lf. = u{'^h,{l + l5Hiy') (B.2) 

Here, we keep the incoming particle is always a massless neutrino, while the outgoing 
can be either a neutrino or a massive lepton. For anti-lepton reaction, 

l^ = v(-u')j^(l + ^,)v{-u) (B.3) 

For completeness, the needed matrix element for /5~ decay is: 

l, = u{eh^{l+^,)v{-u) (B.4) 

The calculation of the cross-section include the expressions /a/*. For the evaluation 
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of this expression, the following facts concerning the 7 matrices are helpful, 

{7^7^} = 2^'^'^ 
{75,7^ = 
(1 + 75)' = 2(1 + 75) 
707^70 = lu 
^?^{7A7A«7i^7M'} = 4:{9x^^9u^,' - gxugi^f,' + 5'A/.'5'/.i.) 
^^{7a757m7^7m'} = ~^'i'^\^lu^l' 



where the metric tensor g 



flU 



/ 1 \ 
0-100 
0-10 



\ -1 / 
For any two spinors tpi and 1^2 and any 4x4 matrix F, 



{iJiTtp2y = V^27or7o^i 



thus, for lepton reaction, 



hi: = U^'''\k2h,il + 75)w^"^H^l)w^"^H^l)(l - l5huU^'''\k2 



(B.5) 



To obtain the cross section it is necessary to sum over initial and average over final 
lepton helicities. 

^ = ^ E ^^"'^^2)7^(1 + 75)w^"^^(A;i)m("^)(A;i)(1 - 75)7.w^"'n^2) (B.6) 

aia2 



It is helpful to use the fact that. 



J2nf{k)uf{k) 



kfj.'y^ + '^ 



2K 



P /gg' 



(B.7) 
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Thus (the mass term falls due to 7 matrices algebra), 



z/^iy'^' 



finally, 



Ut = ^^tr{7A(l + 75)7m7.7m'} (B-S) 



hit = ^{^^< + ^'x^'^^i^^'^ - i^xi^u^'^^^'^'} (B.9) 



the calculation for anti-neutrino reaction is similar and gives 

Wu = -^i^^^'u + ^x''u'y^.^'^. + ^^x^.u^^'l^''l^''''} (B.IO) 

The general neutrino reaction cross-section requires the following expressions, 

TJ* = l + u-i)' (B.ll) 

y* = l-i).i)' + 2{i>-q){u'-q) (B.12) 

y* = q.(0 + ,)') (B.13) 



^{f-f*-hl*,) = l-{u-q){u'-q) (B.14) 



{I X l*)3 = S X 2iq ■ {u - i^') (B.15) 

where g'^ = z/'^' — z/'^ is the momentum transfer. S* is +1 for neutrino scattering and 
l3~^ decay, and —1 for anti-neutrino scattering and /3~ decay. 



NON-RELATIVSTIC EXPANSIONS 
OF DiRAC SPINORS EXPRESSIONS 



Dirac spinor of a mass M, takes the form: 






where Ep = ^p^ + M^ is the energy of the particle. It is useful to recall that u = u^'^q. 
This convention leads to a normalized density, i.e. u^u = 1 
In this representation, 

where we use 7^757'' = —75 and 7^757 = — (x. 

One can now expand the needed matrix elements in the inverse mass: 

u{p', a'h%{p, a) = x'l{l-^,+ ^q■ ^^) X. (C.3) 

u{p\ (t')^u{p, a) = ^Xa{P -iff^ q)Xa (C.4) 

u{p', a')757^(P, cr) = -^x'K^ " P)Xa (C.5) 

m(P , (^'hMP^ ^) = '^Xai^ ■ Q)Xa (C.6) 

HP , (T h^lHp, a) = -xl K 1 - ^777? + ^-rvJT + 7775 I Xa 



M2 / 4M2 8M2 

(C.7) 



where P = p + p' = — i(V — V) and q = p — '^ . 
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Feynman Diagrams for xPT 



The Feynman rules are exhausted from the interaction Lagrangian, and from the 
currents. This is done by expanding the complete term to first order in the pion field. 
The diagrams are then easily achieved. 

It is valuable to remind here the pion propogator p^^, and the nucleon pro- 



pogator 



y-M- 



D.l Leading Order Diagrams 

The pion-nucleon interaction lagrangian to leading order in the pion field is given in 
Eq. (3.17). The resulting Feynman rules appear in Fig. D.l. 



jf-To Ar^ 



q 



p 



1 ^ahc /,^c 



^abc ^^c 



k" 



Figure D.l: The lowest order pion - nucleon Feynman diagrams. Nucleons are 
indicated by solid lines, whereas pions are indicated by dashed lines. 
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i_ oif ,,t: 



e""7' 



■tgAj'-rsY 



f^k"'uS''' 



ail 



aiJ,0_ 



Figure D.2: Vertices in the axial current. Crossed circles indicate the attachement 
to the external probe. 



The axial current to the same order is given in Eqs. (3.30 - 3.31), and appear in 
Fig. D.2. 



To lowest order in the pion field, one can get the following vector vertices from 
Eqs. (3.27 - 3.28): 



D.2 Next-to Leading Order Diagrams 

The NLO pion-nucleon interaction lagrangian to leading order in the pion field is given 
in Eq. (3.21). The resulting Feynman diagram is identical to the diagram drawn in 
the right part of Fig. D.l. The strength of the vertex is: 



2i 



{d,e-''jKq.cr^^'' - csk^q'-d'^''] 



(D.l) 



The NLO contribution to currents within the nuclei appear in Eqs. (3.31 - 3.28). The 
resulting vector current includes two pions, thus of higher order. The resulting axial 
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Figure D.3: Contact axial vertex in NLO. 

vector diagram is identical to the one which appears in the center of Fig. D.2. The 
additional strength is: 

{-c^e'^'^^qucr'"' + CsQ^S'''] (D.2) 



The NLO Lagrangian includes in addition contact terms, as given in Eq. (3.22). 
The form of the vertex is a contact term with an outgoing pion line, and the strength 
of the resulting vertex is —^^^75 /iTa- 

The contact term contribution to the vector currents is zero in the non-relativistic 
limit. However, the contribution to the axial current diagram does not vanish Di'j^'j^y 
with the diagram in Fig. D.3. 



D.3 Example: Single Nucleon Operators For Neu- 
trino Scattering 

One can use the Feynman diagrams of the previous section to draw neutrino scattering 
on a single nucleon: 
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afi. 



® 



fc'' 



-® 



ail 



Figure D.4: Leading order 1-body axial vector current. 

The resulting strength is identical to that achieved by standard nuclear physics ap- 
proach (SNFA): 



M'^^^{l)=igAu{p')^Ay 



(D.3) 
The second equality uses the known SNFA axial form factor Fa = —Qa-, and pion 
form factor Fp = —j2i^- By non-relativistic expansion of the nucleon spinor, as ex- 
plained in Appendix C, one gets the single nucleon scattering operators of Eqs. (3.43 - 
3.46). 



D.4 Example: Axial MEC Operators For Neutrino 
Scattering 

In order to continue to axial currents that involve two nucleons with a pion exchange, 
we will find the amplitude for pion production. This is done in Ananyan et al. (2002), 
and the resulting Feynman diagrams are given in Figure D.5. This is however wrong, 
as the first three diagrams (starting upper left side) include pion exchange between 
nucleons without any interaction of this pion with the external probe. The correct 
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/q" 



X 



"1^ 



®' 



yql' 



®' 



® "'' 



W" 



. ql' 
V kl' 



Figure D.5: Pion-production on a nucleon. 



pion production amplitude is thus, 



r75 /i^i 



j^fi ^v 



b^a 



2 



(/+ (^-M) " ' P-m2 2 

27^^ U(p) 



A;^ — ml 



(D.4) 



The two nucleon axial current is achieved through the exchange of a pion, as observed 
in Figure D.6. The amplitude is: 



M'^^(2) = M^^^in 



7 (1 A T 

^—^ n~r^{'P'i) ilb — u(p2) + momenta permutations (D.5) 



Taking the non-relativstic limit (0(1/M)) and the soft pion limit (g, k — > 0), we get 
for the pion production, 



2/^M'^'"^(7r) ^ -m(p')e"'''rV^(p) 



(D.6) 



hence, 



M^ 



-f^(r« X r(^))«(zA + a, x g + a, x h)^^ + (1 - 2), (D.7) 

o^^^/,; ki + ml 
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afi(^M{nJ 



-q^ 



In 

pi 



Pi 



Figure D.6: Pion exchange contribution to the scattering amphtude 



M' 









:i^2) 



(D.8) 



where Pi = p + ff, ki = pi — p*, and g'is the neutrino's momentum. 



An additional contribution to the pion production amphtude comes from NLO (cal- 
culated using the diagrams in Figure D.7): 






2c, 



UM 



^_^a,. I g. 



kf'k 



k^ 



qMP>^u{p) + 
mi 2 



+ i^e"^^ <) q^ 



kf" 



k-q} 6''%{p')u{p) 



(D.9) 



using Eq. (D.5) we can now get the NLO contribution to the pion exchange currents: 

k^'kx 






9x 



/j2 _ j^2 



quQa 



T 



2 g^ — m^ 2 



2C3 

M 



k'' 



q' 



k"^ — m: 



-k ■ q 



2 1^2 («(Pi)«(Pi)«(P2)7''75ym(P2) ) [>(D.10) 
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a/i 



.hq^ 



/ 



aji 



1" 



/ 



< 



hq^' 



Figure D.7: Pion production contribution from NLO lagrangian. 



The non-relativistic reduction leads to no contribtion to the axial charge density and 
to the following expression for the axial current: 



M^Lo(2) 



^9a 



2M/2 
If we combine this with Eq. (D.7), we get 



h (r(i) X f(2)) ai X ^2 + 2c^f^^H, 



^2 ■ ^2 

H + ml 



(D.ll) 



+ Q+C4VfWxf(2))Vix4 + 

+i(.(^)x.^(^)mx,-)}^ + (i^2), 



(D.12) 



The contribution of the contact Lagrangian vanishes for the axial charge operator in 
the non-relativistic limit. Its contribution to the MEC reads: 



M^ ^ iDi ai 



-(l)a 



r 



■(2)a 



Cr2- 



(D.13) 



In order to get the scattering operators, one has to transform to configuration space. 
After doing thus, one gets the final form of the MEC operators, appearing in Eq. (3.48). 




Reduced Matrix Elements 



OF DIFFERENT OPEARATORS 



In this chapter, I give a hst of the needed matrix elements when calculating in the 
HH expansion, as implemented in the program. In the following S is a spin spherical 
tensor operator of order 1. 

MEl 

{K{ls)j\\jjiqx)Yj{x)\\K'{l's')f) = 

^ I J I' 



V(-)'^^'^'^V(2j + l)(2j' + l)j [ 



X 
s J 



{Kl\\jj{qx)Yjix)\\Kr) (E.l) 

using the fact that. 



.1 /(2/ + l)(2J+l)(2/ + l) I J V 



{i\\yAmi') = {-)\ r " 4, " ' M Q ' ^"^-'^ 



We finally get, 

{K{ls)3\\3j{qx)Yj{x)\\K'{l's')3') 



SsA-y^''""'^ 



{21 + l)(2j + 1)(2J + 1)(2/' + l)(2j' + 1) 



X 



47r 

[ ^ ^'\ili [: {ki\jjM\k'i') (E.3) 
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ME2 

We use the fact that S ■ Yjlm{x) = [S^^) ® F^jj^^ 

{K{ls)j\\jj{qx)^ ■ YjLMmK'{l's')f) = V(2j + 1)(2J + l)(2j' + 1) 

{/' s' f 
L 1 J }{Kl\\jj{qx)YL{x)\\K'l'){s\\^\s') (E.4) 

/ s J 

putting now Eq. (E.2), we get 



{Kils)j\\jj{qx)i: ■ YjLMmK'{l's')f) 




' (2j + 1)(2/ + 1)(2J + 1)(2L + 1)(2/ + 1){21' + 1) 

47r 

X I ^ ^ ^ ) {sm\s'){Kl\\jjiqx)\\K'l') (E.5) 



ME3 



{KI\\Yjl{x)-W\\K'1') 



.+1, /(2/ + 1)(2J + 1)(2L + 1)(2/' + i: 



I L /' - 1 ' 



X 








/ L I' + l 



4 ' ^' ^\y ' ° /i±i(;,,|M-n|,,r)] (E.6) 

1 /' + 1 L 1 I //'+l 1 /'\ 2/' + l^ ^\dr r)^ ^^ ^ > 



I'+l 1 /' 
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ME4 



{K{ls)j\\YjL{x)-V\\K'il's')f) = 

= 5.,.K-)'+^'+'+V(2j + l)(2j' + l)| [ ^ ^' ){Kl\\Yj,{x)-V\\K'l') 



5...A-V- 



+,+,+1 /(2j + l)(2/ + l)(2J+l)(2L + l)(2/ + l)(2/' + l) I / J V 



Att 



J S J 



I L I'-l 



\ I'-l L 1 j I' 1 /'_i\ 2/' + l^ ^\dr r J^ 

\ y 



/ L I' + l 



+ 



\ I' + l L 1 J I' + l 1 r\ 2/'+l^ ^\dr r)\ ^^ 
\ 0/ 



(E.7) 



ME5 

We note that 



YlkAA)^ ■ S = 5^( 



\L-J 



2J+1 
2L + 1 



Yl®V 



(E.8) 



J M 



Thus 



{K{ls)j\\Y,{x)V ■ ^\\K'{l's')j') = ^(-)^-V(2j + l)(2J+l)(2/ + l) 



/' s' f 

x{s\\t\\s'){ J 1 L }{Kl\\YjLix)-\/\\K'l') 
I s j 



(E.9) 
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Radial Integrals 

The radial matrix element to be calculated is, 

{{KN_,lN)Kr,\d{LqpsmeN)\{KN_,l'j,)K'^) = 

■K 

Jo 
X sin''^ On cos^^-^ ^7vi'i"^Hcos20Ar) 

xd(Lgpsin^7v) sin'^ On cos^^-^ ^^Pif'^'^(cos2^7v) (E.IO) 

where N = A-l,n= ^'n-^^-^-^'n ^ a = /^ + i, and /3 = K^v-i + ^• 

/o = {{KA-ilN)KN\jj{LqpsmeN)\{KA-A)K) (E.ll) 

h = {{KA-ilN)KM\jj{LqpsmeM)- ^-—\iKA-A)K) (E.12) 

Lqp sm Un 

h = {{KA-ilN)KM\jALqpsm9r,)-jr-^-^^\{KA-il'N)K) (E.13) 

a[LqpsinUN) 

using the fact that 

— ^ (sin'^ ^;v cos^^-^ ^^Pi"^)(cos20^)) = 

= sin''^ ^^ cos-^^~-^ ^7v^i°'^ncos2^Ar) 

/at Kat-i sin6'7v 2n(n + /3) „ sin 6^ at 

sin 9]y cos^ 6^^? (2ra + a + /5) sin 9^ cos^ 6'Ar cos^ 9^^ 

_ 2in + a){n + P) ^.^,^_, ^^ ^^^,^_^_, ^^pM)(^^^2^^) (E.14) 

2n + a + p 

We can now calculate numerically the integrals above. 

E.l 2-Body Matrix Elements 

The 2-body matrix elements needed for the evaluation of the axial MEC include 
several operators, whose reduced matrix elements are not completely trivial. 

A much needed equalities are the following reduced matrix elements, important 



Appendix E. Reduced Matrix Elements of different opearators 118 



also for the isospin part, of pauli matrices: 

{{\l)S'\\a^ + a,\\{ll)S) = Ss,s'VSiS+l){2S + l) (E.15) 

((^^)0||ai - ^2||(^^)1) = -((^^)l|l'?i - ^2||(^^)0) = -2V3 (E.16) 

{(IIm^i X ^2||(^^)1) = ((^^)lPi X ^2||(^^)0) = -2V3^. (E.17) 

These are used both for the spin and isospin operators. 

By making use of the known spherical harmonics coupling: 



H-l,jMh] h k I 



[y(M ^ y(.)]« = i-y^-^^^l^ - - ^ ) r(o (E.18) 



(here [/] = 2/ + 1) we can achieve the following formula, 

{{L'S')J'\m ® {a^ a,Y' Y^YL\\iLS)J) = 
. ... ...^../\T/]\.T']\T.]\T][J] 

/ T" T T \ ( T" T T 

Eli" 



47r 

L" L ^ ^ J ^" L L ) f L" 1 L' \ ( L" 1 V 

j 1 J 5' J' f \ j I S' J' S 



Thus, 



6% ■ YjLM = (r(i) f(2))-[SQ r^li^J) (E.19) 

f^"-r..M = (r«0r(^)r|(-)^V4;^f J I 'l) [^o^Y.rYjM- life ^yS^ 

(E.20) 

where Sq = Si S2, for = x , +, — . 
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